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INTRODUCTION 


The  need  to  investigate  a  large  number  of  free-stream  and  geometric  combinations  for  preliminary 
design  studies  requires  rapid/accurate  computational  models.  It  is  assumed  that  the  flow  field  computation 
can  be  separated  into  inviscid,  boundary  layer,  crossflow  separation,  and  separated  base  region  comput¬ 
ations. 

The  bulk  of  the  theoretical  effort  was  devoted  to  developing  an  inviscid  computational  model  for 
configurations  with  nonaxisymmetric  bodies  that  would  be  valid  at  low  angle  of  attack.  The  objective  was  to 
develop  a  code  for  the  free-stream  Mach  number  range  of  1.2  S  3.0  with  an  accuracy  approaching  that 
of  the  nonlinear  Euler  or  potential  models,  but  which  requires  significantly  shorter  run  times. 

Linear  surface  singularity  methods  l>2  have  been  highly  developed  for  complex  configurations, 
including  high  angle-of-attack  vortex  modelling.  Set  up  and  run  times  are  fairly  long.  Nonlinear 
compressibility  is  not  well  modelled. 

For  axisymmetric  bodies  at  low  angle  of  attack,  component-buildup  methods  can  be  utilized  for  a 
restricted  range  of  configurations  and  free-stream  conditions  that  are  applicable  to  many  current  designs.3-4 
For  nonaxisymmetric  bodies  with  inlets  and  more  general  free-stream  conditions,  the  component  buildup 
approach  is  less  applicable. 

The  Van  Dyke  second-order  potential  method  5  was  adapted  for  the  body-alone  inviscid  prediction 
described  in  References  3  and  4.  The  method  only  requires  the  evaluation  of  axial  source-sink  and  doublet 
distributions  that  are  evaluated  as  the  solution  is  marched  down  the  body.  Extension  of  the  model  to 
nonaxisymmetric  bodies  requires  evaluation  of  the  first-order  solution  in  the  flow  field.  The  axisymmetric 
body  solution  at  angle  of  attack  is  not  a  full  second-order  solution. 

The  second-order  Van  Dyke  model  has  been  extensively  developed  for  nonaxisymmetric  bodies  with 
planar  discontinuities  including  inlets.  ®-8  In  this  report,  extensions  to  bodies  with  blunted  noses  and  lifting 
surfaces  are  given.  The  numerical  techniques  are  finite-difference  implicit  and  explicit  marching.  The 
current  method  tracks  the  origin  of  disturbances  or  nose  Mach  cone.  This  limits  the  small  disturbance 
solution  to  bodies  that  lie  within  the  Mach  cone.  The  work  described  here  and  in  References  6-8  required 
original  research  that  did  not  exist  in  the  current  literature. 

At  higher  Mach  numbers,  local  pressure  solution  methods  9  are  utilized  (generally  for  Mach  numbers 
greater  than  3.0).  However,  for  better  accuracy,  an  Euler  solution  such  as  given  by  the  SWINT  code  l0-11  is 
recommended. 

The  effects  of  viscosity  (e.g.,  surface  skin  friction,  crossflow  separation  at  high  angles  of  attack,  and 
base  wake  flow)  are  treated  by  rapid,  very  approximate  methods. 
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2.0  GEOMETRY,  FREE-STREAM  VELOCITY,  AND 
FORCE  CONVENTIONS 


2.1  BODY  GEOMETRY 

The  body  is  assumed  to  have  a  vertical  plane  of  symmetry.  For  the  linearized  computations,  the  body 
(including  fins)  must  lie  within  the  Mach  cone  of  disturbances.  For  a  blunted  body,  the  origin  of  the  Mach 
cone  begins  upstream  of  the  nose.  However,  body  geometry  is  always  referred  to  the  nose. 

A  half-plane  for  a  general  body  is  depicted  in  Figure  1.  Body  geometry  is  given  in  cylindrical 
coordinates  and  the  body  radius,  rt,,  must  be  single  valued  in  the  azimuth  angle,  9.  For  the  blunted  body 
shown,  the  choice  of  origin  of  coordinates  0’  will  be  discussed  in  Section  3.7.  The  body  is  assumed  to  be 
continuous  between  planar  discontinuities.  However,  rb  is  assumed  to  be  continuous  in  0. 

An  initial  blunted  section,  0  to  the  plane  00,  must  be  circular  in  cross  section.  The  rb  variation  with  x 
need  not  be  that  of  a  circular  arc.  Typical  planar  discontinuities  are  shown  for  00,  AA,  and  BB:  00  may 
have  a  slope  discontinuity  or  a  radius  of  curvature  discontinuity;  AA  depicts  a  radius  or  planar  inlet 
discontinuity;  and  BB  depicts  a  slope  discontinuity. 

The  body  is  divided  into  sections  equal  to  the  number  of  planar  discontinuities  plus  one.  The  body  is 
assumed  to  be  piecewise  continuous  to  second  order  in  x  derivatives  between  planes  of  discontinuity.  No  limit 
to  the  number  of  sections  is  assumed.  Any  section,  except  the  first,  may  have  a  set  of  fins  (up  to  six  fins  per 
set).  Each  section  may  be  similar  or  nonsimilar.  Similarity  implies  a  basic  cross-sectional  outline  that  is 
scaled.  An  axisymmetric  body  is  similar  throughout. 

Each  cross  section  is  quarter-  or  half-plane  symmetric  and  dr^/dQ  is  assumed  to  be  zero  at  the  ends  of 
the  symmetry  intervals. 

Body  geometry  is  input  analytically  by  code  source  updates  or  numerically  as  planes  of  rb  vs  0  point 
sets.  Each  section  may  have  a  different  mode  of  body  geometry  description. 

For  similar  sections,  only  one  reference  x  plane  of  rb,  3rb/30,  and  d2rb/302  need  be  evaluated  for  the 
quarter-  or  half-plane.  For  the  remaining  x  locations,  rb,  dryjdx,  and  32rb/3x2  are  required  only  at  the  initial 
computational  0  plane.  For  the  numerical  mode  of  input,  rb,  3rb/3x,  d2rb/3x2  are  computed  using  a  four-point 
Lagrange  polynominal  fit;  drb/30, 32r/302,  and  32rb/3x30  are  computed  numerically  using  third-order  splir' 
fits. 


2.2  FIN  GEOMETRY 

The  fins  are  treated  as  thin  fins  in  a  manner  very  close  to  that  described  in  Reference  10.  The  fins  are 
assumed  to  be  symmetric  about  a  midplane  that  lies  on  a  constant  0  plane.  Figure  2  depicts  the  sectional 
view,  CC,  from  Figure  1.  Also  depicted  is  an  intermediate  computational  plane  and  an  expanded  view  of  the 
fin  cross  section. 

Again,  symmetry  is  built  into  the  code.  The  fin  configuration  is  input  as  quarter-  or  half-plane 
symmetric.  Thus,  if  the  body  has  quarter-plane  symmetry  and  the  fin  set  has  half-plane  symmetry,  the  body 
information  is  for  O<0<n/2  and  the  fin  information  is  for  two  fins  as  shown  in  Figure  2. 
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As  indicated  in  Section  2,1,  up  to  six  fins  per  full-plane  are  allowed.  For  fin  quarter-plane  symmetry, 
only  one  or  two  fins  are  allowed  in  the  interval  0  <Q  <  n/2.  Two  fins  are  allowed  if  one  of  the  fins  lies  on  either 
0  =  0  or  0  =  n/2  planes. 

For  fin  half-plane  symmetry,  one  to  four  fins  arc  allowed  in  the  interval  0<9<n;  four  fins  are  allowed 
if  two  of  the  fins  lie  on  0  =  0  and  0  =  n  planes.  None  of  the  fins  for  fin  quarter-  or  half-plane  symmetry  need  to 
be  identical. 

The  indices  i,  j,  and  k  refer  to  a  computational  grid  that  is  not  shown  in  detail  in  Figure  2  and  will  be 
discussed  later  in  this  document.  £,  is  a  shearing  transformation 


_  r-rb  (2.1) 

that  maps  the  region  between  the  body  and  the  Mach  cone  to  a  rectangular  region. 

The  thin-fin  assumption  implies  that  the  fin  surface  lies  close  to  the  fin  midplane.  For  computational 
purposes,  the  fin  is  assumed  to  have  sharp  edges.  Fin  boundary  conditions  are  evaluated  at  the  midplane. 
Thus,  points  A  and  B  coalesce  to  two  computational  points  on  either  side  of  the  midplane.  The  fin  plane  is 
located  by  the  index  j=JFIN.  The  9f  values  are  obtained  from  the  body  input  array;  8  j ,  02, .  •  9 {])  ■  • ,  9m  = 
0(JM);  to  be  discussed  later. 

From  the  geometry  of  Figure  2, 

(9-9} 

qx  =  0  -  Qf(JFIN)  -  o  v  1 


o  =  tan  1  (Mr), 


(2-3) 


where  tp,  o,  and  fare  assumed  to  be  small. 

All  moveable  control  deflections  are  treated  only  as  streamwise  slopes.  For  a  fin  deflection,  in  reality, 
leading  and  trailing  edges  are  displaced  out  of  the  original  fin  plane.  The  location  of  the  deflected  fin  hinge 
line  is  not  an  input.  Individual  fin  deflections  are  input.  A  positive  deflection  occurs  when  the  trailing  edge 
of  the  fin  lies  clockwise  from  the  fin  plane. 

A  consideration  offree-stream  conditions,  fin  symmetry,  fin  deflections,  and  body  symmetry 
determines  the  values  of  0i  and  0max  and  whether  the  computation  is  quarter-,  half-,  or  full-plane.  This  will 
be  discussed  later. 

Figure  3  depicts  the  geometry  used  to  describe  a  single-fin  planform.  XMI  is  the  minimum  value  of  x 
for  the  fin.  XMA  is  the  maximum  x  distance  of  the  fin  from  the  origin.  RMI  is  the  minimum  radial  distance 
from  the  x  axis.  s  +  RMI  locates  the  general  radial  distance  of  a  fin  point.  SMA  +  RMI  is  the  maximum  radial 
extent  of  the  fin.  All  the  values  in  this  paragraph  are  input  numerically.  With  this  type  of  input,  the  same 
fin  may  be  positioned  any  place  on  the  body. 

The  planform  outline,  its  inclinations,  and  the  fin  thickness  derivatives  are  given  as  analytic  functions 
of  s  and  x.  The  planform  outline  is  given  as  inner,  SI,  and  outer,  SO,  values  of  s.  Currently,  the  fin  surface  is 
assumed  to  be  a  continuous  surface  for  most  applications.  For  the  high  Mach  number  range,  discontinuous 
fin  slopes  may  be  input  as  is.  For  the  potential  flow  computation,  discontinuous  slopes  are  scaled  down;  fin 
axial  force  is  then  scaled  back  up.  The  scaling  is  also  needed  for  thick  continuous  surfaces. 
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Also  shown  in  Figure  3  are  the  traces  of  constant,  L,  values  on  the  I'm.  For  the  high  Mach  number 
regime,  these  traces  are  based  on  a  pseudo  Mach  cone  that  lies  outside  all  fins  in  the  section.  The  edge  where 
a  f,  trace  enters  the  fin  is  a  leading  edge.  The  edge  where  a  L,  trace  leaves  the  fin  is  a  trailing  edge.  If  the 
trailing  edge  is  subsonic,  it  is  a  side  edge. 


2.3  FKEE-STREAM  VELOCITY  VECTOR 

The  Cartesian  components  of  the  free-stream  dimensionless  velocity  vector  are  given  by 

qx  —  cosa  cos  (F  i'  +  sin  P '  j'  4-  sin  a  cos  P'£' 


(2-4) 


2.4  FORCE  CONVENTION 

The  axial  force,  F^,  acts  in  the  x  direction;  the  normal  force,  Fn,  acts  in  the  7.  direction;  and  the  side 
force,  Fy,  acts  in  the  negative  y  direction.  The  roll  moment,  M;,  acts  in  the  negative  x  direction;  the  yawing 
moment,  Mn,  acts  in  the  negative  z  direction;  and  the  pitching  moment,  Mm,  acts  in  the  negative  y  direction. 


2.5  SYMMETRY  CONSIDERATIONS  AND  NONZERO  FORCE  COEFFICIENTS 

As  earlier  indicated  in  Section  2.2,  the  range  of  G  and  the  nonzero  force  coefficients  are  determined  by 
the  free-stream  velocity  vector,  body  geometry  symmetry,  fin  symmetry,  and  fin  deflection.  In  the  coding,  the 
individual  force  and  moments  are  labled  1  through  6  for  axial  force,  normal  force,  pitching  moment,  side 
force,  yawing  moment,  and  roll  moment,  respectively.  There  are  four  possible  combinations  or  modes,  LM,  of 
G  range  and  active  force  coefficients.  The  LM  value  may  change  from  body  section  to  body  section  and  is 
associated  with  the  highest  value  of  force  label.  LM  is  an  input  variable.  LM  =  I  is  for  quarter-plane 
symmetry  for  body  and  fin,  a  =  0,  P'  =  0,  and  mirror  image  fin  deflections  about  6=  n/2.  Only  the  axial  force 
is  nonzero.  LM  =  3  is  for  vertical  half-plane  symmetry.  It  is  obtained  bv  various  combinations  where  P'  is 
always  zero.  Active  force  labels  are  1  through  3.  LM  =  5  is  the  upper  half-plane  symmetry  associated  with 
any  combination  of  yaw  deflection  and/or  side  slip  and  body  quarter-plane  symmetry.  Active  force  labels  are 
1 , 4,  and  5  LM  =  6  is  for  the  full  range  of  0  and  forces.  Configuration  symmetry  may  vary  from  section  to 
section.  LM  ^  LML,  where  LML  is  the  previous  section  value  of  LM,  is  required  oy  the  hyperbolic  character 
of  the  flow.  Thus,  LM  may  vary  from  1  to  3,  1  to  5, 1  to  6, 1  to  3  to  6,  1  to  5  to  6,  or  5  to  6  in  the  downstream 
direction. 


3.0  THEORETICAL  DEVELOPMENT  FOR  THE  FIRST-  AND 
SECOND-ORDER  POTENTIAL  SOLVERS 


3.1  FIRST  - AND  SECOND-ORDER  POTENTIAL  EQUATIONS 
The  full  nonlinear  dim^.isionless  potential  equation  is  given  by 

,,2 

1  +  M2  (1  -Q2) 

l  2  00 

A  small  disturbance  from  the  free-stream  dimensionless  velocity  vector  is  given  by 


V-  Q=  ~Q-  V  (Q2) 

LJ 


(3-1) 
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Q  =  < 7,3  +  <7  =  (cos  a  cos  P'  +  u)  i'  +  +  q  '3' 1 

where  Q  is  the  total  velocity  vector,  q  is  the  disturbance  velocity,  qa  is  the  free-stream  velocity,  and  q  and  q 
are  crossflow  components  of  q^  and  q,  respectively. 

Substituting  Equation  (3-2)  into  (3-1)  and  neglecting  higher-order  terms,  one  obtains  the  well-known 
first-order  wave  equation 


V  • 

c 


lc 


dll, 

p2  — -  =  0 

ax 


(3-3) 


Vc  is  the  crossflow  gradient  operator.  The  subscript  i  indicates  first  order.  Equation  (3-3)  is  supplemented  by 
the  two  crossflow  irrotationality  equations. 


(VXQl)c  =  0 


(3-4) 


An  improved  solution  to  the  first-order  problem  is  obtained  by  evaluating  the  neglected  nonlinear 
terms  using  the  first-order  solution  and  solving  the  nonhomogeneous  wave  equation  for  a  second-order 
iterative  solution.  Second-order  equations  are 


V  ■ 

e 


l2 c 


Q 

c 


M 


-  2-y  2  v 

—  -  a,  V  •  q 

2  g  4  lc  "  4 


lc 


(3-5) 


(Vxq2)c  =  o 


(3-6) 


a=ui  <*«?- 


+  cos  a  cos  P' 


2  ,  Y  “ 1  „  2  2  , 

q,  +  - M  u,  +  q  •  q 

M  9  oo  1  ’1 


(3-7) 


n 


7lc +  (9?  + ?i- 


q 

■  oo  * oop 


+ 


P2 


+  qlc  x  (qu 


(3-8) 


The  subscript  2  refers  to  the  second  order.  The  first-  and  second-order  equations  are  in  nearly 
conservative  general  vector  form.  Only  the  last  term  on  the  right  hand  side  of  Equation  (3-5)  is  in  non¬ 
conservative  form. 

The  boundary  condition  on  all  surfaces  is 


Q-  n  =  0 


(3-9) 
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where  n  is  the  outward  normal. 

Any  set  of  coordinates  may  be  used  for  the  crossflow  plane.  An  ideal  set  of  body  fitted  orthogonal 
coordinates  would  transform  the  region  between  the  body,  including  fins  and  the  bounding  Mach  cone,  to  a 
rectangular  computational  region.  This  would  require  a  numerical  solution  of  Laplace's  equation.  This 
would  lead  to  a  simplified  boundary  condition  evaluation.  However,  the  disturbance  envelope  Mach  cone  is 
circular  and,  hence,  the  choice  of  cylindrical  coordinates  eliminates  the  grid  generation  problem.  However, 
the  range  of  configurations  is  limited  unless  one  introduces  the  complex  multizone  approach  of  Reference  12 


3.2  COMPUTATIONAL  COORDINATES  AND  GRIDS 

The  shearing  transformation  has  been  introduced  in  Equation  (2-1).  The  basic  variables  £,  and  9  are 
further  clustered  independently  for  unequally  spaced  grids. 

The  9  grid  may  be  uniform  or  nonuniform;  9  points  are  clustered  near  points  of  greatest  curvature  of 
the  body.  It  is  assumed  that  the  nonuniform  grid  9  is  mapped  to  the  uniform  unspecified  grid  4>.  For  a 
uniform  grid,  only  the  maximum  9, 9m,  compatible  with  the  body  symmetry  is  input  and  divided  equally 
according  to  a  maximum  input,  j  —  JM.  For  a  nonuniform  grid,  the  individual  9  computational  points  are 
input.  Transformed  first  and  second  derivatives  are  evaluated  by  applying  central  or  one-sided  differences  in 
Combinations  that  are  needed  are 


a/30(  )  = 


)  ~ 


l 


A<t> 


392 


( 


l__S 9 

dtp2  dty1 


de 

*** 


(  ) 


1 


d29 

d§2 


2  dQ 

(A$)2/(—  A4») 
dty 


(3-19) 


(3-11) 


The  quantities  d9/d<J>  A4>  and  d29/d(}>2  (A<J>)2  are  evaluated  by  central  or  one-sided  difference  approximations  of 
O  [(A<}>)4]  accuracy.  Ct  and  C2  are  related  to  the  finite-difference  approximations. 


The  £grid  may  be  uniform  as  determined  by  the  input  KM  value  or  nonuniform;  £  1  =  £(1)  =  9.  The 
nonuniform  grid  may  vary  as  a  geometric  progression  in  which  £,  2  =  £  (2)  =  DO  or  may  be  input 
individually,  ^clustering  is  necessary  far  downstream  on  long  bodies  or  for  small  span  fins.  The  geometric 
progression  clustering  grid  is  given  as 


^  =  DO 


L=  k-2 

V 

/_ 

L~0 


L  _ 
a  = 


DO  ( 1-g*-1 ) 
1  —a 


k  =  2,3 . KM+ 1 


a  is  determined  at  k  =  KM  +  1  from 


(3-12) 
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,  =  OOd-g™)  (3.,; 

1  -  a 

a  is  determined  by  Newton-Raphson  iteration;  a  &  1.  Once  a  is  determined,  Equation  (3-12)  evaluates  the 
geometric  progression  distribution.  A  transformation  from  5  to  an  unspecified  uniform  5  grid  is  applied  in  a 
manner  similar  to  that  for  the  9  gridding.  First-  and  second-order  derivatives  similar  to  that  in  Equations 
(3-10)  and  (3-11)  need  to  be  evaluated.  Central  and  one-sided  difference  approximations  of  0  [(V<)4]  are  used. 


3.3  VELOCITY  ADVANCEMENT  EQUATIONS  IN  5,0  COORDINATES 

The  first-  and  second-order  equations  may  be  recast  into  an  almost  conservation  form  using  5  and  0 
coordinates. 

The  first-order  equations  are 


a/05 


0/0X 


R  (-  p2  Ex  u  t  -I-  ot)  +  Et  wx 

+  a/ax 


+  a/ae 


—  p2x2/? (l/p  -Rb)u{ 


x(l/p  —  Rft)  wl 
=  0 


x(l/p  - 


+  a/a5 


—  u,  +  E  v 
l  z  l 


=  0 


(3-14) 


(3-15) 


a/ax 


x2R  (1/p  -  Rb)u>l 


+  a/a5 


RE  w  -  E  u 

x  1  T  1 


-a/ae 


xd/p  -  RJ  Uj 


=  o 


Here,  R  =  r/x,  Rb  =  rb/x,  and 


(3-16) 


E 

X 


dR. 

R  +x—  (1-5) 

0X 


Er=  -(1-5) 


uj,  vlf,  and  wj  are  cylindrical  coordinate  velocities. 


7 


NSWC  TR  86-253 


The  second-order  equations  are  of  similar  form. 


d/dx 


-pVfl(l/p-R6)u2e 


+  a/ae 


*d/p  -u6)% 


+  d/a(, \x 


«<~P  ExU2e  +  V2e)  +  ETW2e 


=  -q 


N 


d/3i. 


x(Rv{  + 


+  a/ae 


xd/p-fl^ 


a/ax 


x  (1/p  — K6)  u2 


+  d/di, 


ExV2~U2 


=  0 


a/ax 


x*R(\/V-Rb)w2 


+  d/di, 


RE  uj0 

x  2 


-£tu2]1 


-  a/ae 


x(l/p-K.)u 


6  2 


=  0 


2  2  (2  — Y) 

Ml 

U2e  =  U2  +  ^ A 

V2e  =  V2-Ml<ir 
W2e  =  W2~Ml% 


A  is  given  by  Equation  (3-7)  with 


qa  =  cos  a  cos  P'  i'-f  ( sin  a  cos  P'  cos  6  +  sinP'  sinQ)j' 

+  ( sin  P'  cos6  —  sin  a  cos  P'  sin  6)  k' 

=  u  i'  +  v  j'  +  w  k' 

co  on  v  on 


i',  j',  k'  are  unit  vectors  for  the  cylindrical  coordinate  system. 


q2  „  (Y-l  )Ml 

q  =o.  —  +  v  (q  •  q.  +  qh  +  - - - 

’r  1  o  oo  ^oo  ^1  m  ft2 


Wl<ll  2  Y  —  1  2 

qa  —  -  +  w  (q  •  q,  +  q,)  +  — ~  M 2 

2  oo"*®  p2  ® 


w 

-~(o21-lv2l)  +  oa>ulw1 


(3-17) 


(3-18) 


(3-19) 


(3-20) 


(3-21) 


(3-22) 


(3-23) 
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The  equations  in  this  section  will  be  utilized  to  advance  the  field  velocities  by  an  explicit  marching  scheme  to 
be  discussed  later. 


3.4  IMPLICIT  MARCHING  SOLUTION 

An  implicit  marching  solution  is  used  for  continuous  pointed  bodies  with  no  lifting  surfaces  and  for  the 
initial  section  of  pointed  bodies.  The  velocity  vectors  qt ,  q2  are  replaced  by  the  gradient  of  a  potential 


<Ie  =  V4><  f  =  1-2 

A  similar  conical  potential  function  is  introduced  as 


(3-24) 


$(  =  xF(  ( x ,  t,,  0) 

The  full  first-  or  second-order  boundary  condition  is  given  as 


(3-25) 


dr.  ,  8Rl 

V,=  —U,+  - - -W. 

{  dr  *  Ru  80  e 

o 


U,  V,  W  are  total  cylindrical  coordinate  velocities. 


(3-26) 


Equations  (3-14)  and  (3-17)  may  be  written  in  nonconservative  form  after  substitution  of  Equations 
(3-24)  and  (3-25)  as 


<W,  +  C2(ft),  +  W,  +  C4(fte),  +  C5«(fu)(  +  2 Ce«(F  ), 
+  C. x\F  )=(RHS),  (=1.2 

0  XX  c  c 

(RHS)l  =  0.0 


(3-27) 


The  boundary  conditions  may  be  written  as 

(Fl)t  =  Bl  +  B2Ft  +  B3^FQ>e  +  B2X^Fx>t 

Equations  (3-27)  and  (3-28)  could  be  written  in  £,  variables  as  well. 


At  the  Mach  cone,  the  conical  potential  functions  are  known  to  be  singular  and  assumed  to  behave  as 


F,  =  (l-9' 


3/2 


Sll  +  S12^  W* 


(3-29) 


F2  =  s21a-u  +  s22(i-y3/2 


(3-30) 


Second-order  accurate  difference  schemes  are  the  goal.  For  k  =  2,  KM- 1,  central  difference 
approximations  are  utilized  for  all  C,  derivatives.  Central  differences  are  utilized  for  all  <|>  derivatives  as  well. 
At  symmetry  planes,  Fe  and  Fyj  are  set  to  zero.  At  k  =  KM,  £  derivatives  are  obtained  by  differentiating 
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Equations  (3-29)  and  (3-30).  Su,  S12,  S21.  and  S22  are  computed  by  evaluating  Equations  (3-29)  and  (3-30)  at 
k  =  KMandk=KM-l. 

Only  one  previous  plane  of  data  is  needed  in  evaluating  the  x  derivatives.  From  a  Taylor’s  series 
expansion  about  a  known  plane  of  data, 


2 

It  is  assumed  that  Fxx  is  constant  for  the  interval  Ax.  Here,  Ax'  =  x  —  x;. 


F<  =  F<,i  +  <Fx\*x'  + 


(Ax)" 


F  -  F 
t,i  + 1  f ,  i 


2  (F  ).  . 

X  C,  l 

Ax 


Differentiating  Equation  (3-31)  and  substituting  Equation  (3-32)  yields 


(3-31) 


(3-32) 


(*V#,  +  l=-<FA,+  £ 


F  -  F 
t,i+ 1  e,i 


(3-33) 


(Fx)f  and  Ff  are  stored  at  i,  xj.  The  solution  is  desired  at  i  +  1,  x,  + Ax.  These  x  derivatives  then  maybe 
substituted  into  Equations  (3-27)  and  (3-28).  At  the  initial  plane,  (F*)^  1  and  (Fxx)f,  1  are  set  to  zero.  At  x  =  0, 
the  flow  is  conically  similar.  The  advantage  of  an  implicit  scheme  is  that  the  initial  plane  of  data  is 
established  immediately  and  large  steps,  Ax,  may  be  taken. 


At  k=  1,  the  potential  equations  are  evaluated  with  the  aid  of  the  boundary  equations.  F^x  is  obtained 
by  differentiating  Equation  (3-28)  with  respect  to  x.  Initially,  F^e  was  evaluated  by  differentiating  Equation 
(3-28)  with  respect  to  9.  A  more  stable  approximation  proved  to  be 


F^  is  evaluted  as 


<Vu 


F  - F  -F  + 
2J+1  2J-1  IJ  +  1 


l.J-1 


(3-34) 


(F 


tt'U 


2x 

(Ay2 


F  -  F  -(F) 

2  ,j  riJ  ^1/1 J 


(3-35) 


The  resultant  finite-difference  scheme  leads  to  a  block  tridiagonal  matrix  equation  for  (Fi)  j  +  i  and 
(^2)  i+i-  Each  block  is  tridiagonal  for  quarter-  or  half-plane  symmetry.  For  a  full  360-deg  solution, 
additional  wrap-around  terms  are  present  for  each  block.  Solution  of  the  matrix  equations  is  obtained  by 
lowerrupper  decomposition  involving  the  block  matrices  (see  Reference  7  for  details).  Velocities  are  obtained 
from  the  definition  of  the  potential  and  conical  potential  function  combined  with  boundary  condition  and 
finite-difference  approximations  (see  Reference  6  for  details  of  the  velocity  relationships).  Velocities  are 
evaluated  only  at  the  boundary  except  where  a  transition  from  an  implicit  solver  to  the  explicit  solver  is 
required. 
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3.5  PLANAR  BODY  DISCONTINUITY  JUMP  CONDITIONS 

At  the  discontinuity,  locally  two-dimensional  jump  relations  are  sought  that  are  compatible  with  the 
first-  and  second-order  differential  equations.  Since  the  discontinuity  is  planar,  disturbances  are  propagated 
into  the  flow  along  a  surface  that  is  the  envelope  of  the  Mach  cones  emanating  from  the  surface  cut  curve  just 
downstream  of  the  discontinuity. 

Figure  4  depicts  a  cylindrical  coordinate  system  and  the  s,  n  coordinate  system  associated  with  the 
discontinuity  cut  curve;  s  is  along  the  curve  and  n  is  normal  to  it.  The  wave  front  grows  in  the  n  direction. 
Jumps  occur  across  the  normal  to  the  wave  front.  Velocity  components  tangential  to  the  wave  surface  are 
conserved.  The  discontinuity  surface  is  not  tracked  into  the  flow. 

First-order  equations  in  s,  n,  x  coordinates  are  given  as 


1 


dv 


si 


1 


1  +Kn  <5s 


1  4-km 


d/dn 


V  ,  (1  +  Krt) 
ftl 


du, 

-  p2  —  =  0 
dx 


dv 


n  1 

dx 


du. 


dn 


=  0 


d/dx 


(1  +km)  v 

Si 


ds 


(3-36) 


Here,  k,  is  the  cut  curve  curvature.  vnj  and  vsi  are  normal  and  arc  velocity  components.  Jump  relations  are 
established  by  an  application  of  Gauss’s  theorem  for  the  s,  n,  x  system  on  either  side  of  the  as  yet 
undetermined  discontinuity  surface.  This  approach  is  known  as  the  method  of  "weak  solutions’’'^.  Applying 
the  method  of  "weak  solutions”  to  the  set  of  Equation  (3-36)  yields 


9 

Ao  ,  e  +  Av  ,  e  — P  A u,  e  =  0 

si  s  rcl  n  r  lx 

At;  ,  e  -  A u.e  =0 

fi!  x  In 

Av  e  -  A u.e  =0 

St  X  1  $ 

e„,  e*>  8z  are  components  of  the  surface  normal.  Since  A vsj  =  0,  es  =  0. 
From  Equations  (3-37)  and  (3-38), 


(3-37) 

(3-38: 

(3-391 


e  le  =  —  P 

n  x  r 

Aunl  +  PAut  =  0 

Here,  A  (  )  =  (  )  after  -  (  )  before  =  (  )a  -  (  )  b.  Equation  (3-40)  is  the  well-known  result  that  the 
discontinuity  surface  coincides  with  the  characteristic  surface  for  linear  sets  of  equations. 

From  the  geometry  of  Figure  4  and  Avsj  =  0,  one  obtains, 


(3-40 

(3-41 
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Ausl  =  Au>t  cos  0 1  -f  A  u j  sin  0 1  =  0 


A  i>M  =  A  Uj/cosG' 


The  boundary  condition  on  the  downstream  surface  is 


(3-42) 


(3-43) 


where 


+  AOj  =  —  (1-B)u<b+  B(Ulf)+  Auj)  +tan0’(W16  +  Au;1) 


(3-44) 


ton  0'  = 


_L 


Rb  30 


(3-45) 


Equations  (3-41)  through  (3-45)  lead  to  a  solution  for  Auj,  Avj,  Awj.  The  factor  B  allows  for  application 
of  the  modified,  B  =  0,  or  full,  B  =  1,  boundary  conditions.  Van  Dyke 14  showed  that  B  =  1  is  best  for  surfaces 
with  positive  inclination  jumps  in  slope  and  B  =  0  is  best  for  surfaces  inclined  at  large  negative  angles. 

Applying  the  same  approach  to  the  second-order  equation  leads  to 


Au  „e  —  62Au„e  =  p2Au'e  4-  A  v  e 

n2  n  r  2  z  r  x  n  n 


A°„2C*-  AlVn  =  0 


Avs2  =  0 


(3-46) 


(3-47) 


(3-48) 


Au'=  — -  A  A 

P2 


(3-49) 


=  v  , - h  u  (q,  +  q  ■  a.) 

nl  O  no o  M  “ on  ^1 


+  — —  (v-1)  -  v  +  V  ,V  V  , 

-2  12  a°°  si  s®  nl 


(3-50) 


i/4  3 

ao  (2  — y)  2 

—  —  -  —  4*  v  .  v  , 

p2  2  l  3  *x  nl 

Here,  A  v,i  =  0  has  been  used  in  deriving  Equation  (3-50).  Also, 


v  =  v  cos 0'  —  w  sin 0' 

/JOO  CD  00 


v  =  w  cosQ'  +  u  sin  6' 

5OO  (T)  00 


Equations  (3-46)  and  (3-47)  combined  may  be  written  as 


(3-51) 


(3-52) 
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e  /e  =  -B  + 
n  x  r 


(Aun2  +  p  Au0) 


(3-53) 


Au„ 


Ay„2  +  PAu,  = 


(P  A  u'  +  A  o'  e  /e  ) 

n  n  x 

e  /e  -P 


(3-54) 


Ifen/ex  is  set  equal  to  the  first-order  value  of  -P  in  Equation  (3-54),  one  obtains 


Ay«2+PA“2 


(3-55) 


Equation  (3-55)  leads  to  a  set  of  linear  relations  for  Au2,  Av2,  and  Aw2  and  is  compatible  with  Van  Dyke’s5 
jump  relations  for  the  second-order  axisymmetric  case.  An  alternate  approach  is  to  combine  Equations  (3-53) 
and  (3-54),  which  leads  to  a  quadratic  relationship 


Au2  -  Av  0  A o'  -  P2 

n2  n2  n  r 


(Au/ 


+  Au^Au' 


=  0 


(3-56) 


Equations  (3-55)  or  (3-56)  combined  with  Avs2  =  0,  Av„2  =  AV2/COS  0',  and  the  second-order  downstream 
boundary  condition  lead  to  a  solution  for  Au2,  Av2,  and  Aw2-  The  second-order  boundary  condition  is  given  by 


drh 


+  tanQ'  (,W2b  -r  Au>2) 


(3-57) 


Choice  of  sign  for  the  nonlinear  quadratic  solution  depends  on  whether  the  local  turning  angle  is  compressive 
or  expansive. 


The  turning  angle  in  the  n,  x  plane  for  the  second-order  case  is 


DTH  =  tan~l 


cos  O' 


—  tan 


V2bcos0’  -  W2b sin  0’ 


U, 


2  b 


(3-57a) 


Note  that  the  only  difference  in  applying  the  jump  relations  for  surface  slope  discontinuities  and  inlets 
is  that  the  upstream  velocity  vector  at  the  inlet  lip  must  first  be  determined  by  interpolation  in  the  upstream 
velocity  field.  To  continue  marching  downstream  of  the  inlet,  the  entire  velocity  field  between  the  lip  and  the 
Mach  cone  from  the  origin  must  be  determined  by  interpolation. 


Numerical  solution  without  special  damping  or  tracking  of  flow  discontinuities  becomes  difficult 
downstream  of  a  discontinuity.  An  explicit  MacCormack  marching  scheme  was  developed  in  the  hope  that 
the  internal  propagation  of  discontinuities  would  be  captured. 
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3.6  EXPLICIT  MARCHING  SOLUTION  OVER  CONTINUOUS  BODY  SECTIONS 


Equations  (3-14)  through  (3-19)  may  be  used  to  advance  the  field  velocities.  Characteristic- 
compatibility  relations  will  be  combined  here  with  the  boundary  condition'5  to  advance  the  boundary 
velocity  components. 

The  basic  equations  for  combined  first-  and  second-order  may  be  written  as 


u 

X 


+  “*s'" 


uiK+ 


=  R 


Vx+U^x-U^r  =  ° 


,  .  S,  -  -  ucL  =  u 

x  J,  x  r  t  v 


w  +  w. 


=  ue/r 


(3-58) 


(3-59) 

(3-60) 


a  1  dw 

R“  ^2  rp2  ae  p2 


aA  aA 
—  +  E,  — 
dx  x  $ 


+  -  Ld/d  Arq  )  + 

f  r  S  r 


'JV 


-  d/d  0  +  -  ^  d/d  i, 


(3-61) 


■  1 

1  1 

-^ra'\irVl)  + 

-  a/ae  +  -  Ld/dt, 

r  r 

"l 

P' 

Here,  £r>  £x,  and  are  transformation  partial  derivatives. 

Characteristic  eigenvalues  and  eigenvectors  are  sough’,  for  the  x  plane  with  9  as  an  interior  variable. 
The  set  of  Equations  (3-58)  through  (3-60)  are  of  the  form 

(3-62) 


P  + RP  = RH 

X  t 


The  eigenvalues  and  eigenvectors  are  given  by 

L‘ 


R  = 


:  —  x ./ 

i  | 

=  0 

i- 1.2 

1  r 

** 

P2 

'p2^ 

-i, 

i 

0 

r 

X 

-H  ° 

(3-63) 


(3-64) 
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Solving  for  tho  eigenvalues  from  Equation  (3-63), 


(3-65) 


The  choice  of  sign  for  \i  is  associated  with  the  incoming  characteristic  that  is  of  interest.  Corresponding  left 
eigenvectors  are 


L 


l 

l 


(3-66) 


Multiplying  Equation  (3-62)  by  the  left  eigenvectors  and  utilizing  Equation  (3-63),  one  obtains  two 
compatibility  equations. 


(3-67) 


LlP  +  A./Jp.  =  L'  RH 

X  1  (. 


(3-68) 


The  second  compatibility  advancement  equation  of  (3-68)  may  be  simplified  by  utilizing  the  third 
unused  irrotationality  component. 


-  =  Mr  tor)  0-69) 

86 

This  eliminates  the  r  or  ^  dependence.  For  computational  applications,  Equation  (3-68)  is  evaluated  at 
the  surface. 

The  first-order  compatibility  equations  are  of  the  form 


dE..  8G..  5ff,. 

— +  — ii  +  — ii  =C,  i=  1.2 
8x  8^  80  ll 


E..  = 


+  P“t 


11  COS0' 

E12  =  ClVCl=  xRb,COsQ' 

8r . 


’ll-  2  “trv2 

l 


G,,  —  C0Eti,Cq  1  ^  C3  j  C4 
1 


2  P  cosO'  ’ C*  x(l/p-flfc)’~12 


,G„  -  0 


(3-70) 

(3-71) 

(3-72) 

(3-73) 

(3-74) 
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tfn  =  (tonG'  <W<VC5  = 


(3-75) 


H12=  ~Ul~  TU1 


(3-76) 


C12==0-Cll=C6(C7Ul-U'l)  +  C5(C8Ul  +  C3WrC9U,l)'C6 


3tan  0' 


■,c7  =  Psm9', 


Stan  0 
C*=~dQ~ 


'  (  \ 2  sine'  ( 

-  — (ianO  J  ,  cg=  -  ( 


3ton0'  dtanQ' 

Yl  30 ' 

OtJ  W 


Here,  0',  v„j,  and  vni  have  been  defined  earlier  in  connection  with  the  jump  conditions. 
The  second-order  compatibility  equations  are 


^2.  3G2i  ^2,  „  3S  3T 

+  — —  -j-  _  (7  ■+•  a ..  (  +  )  t  — 1,2 

dx  dt,  36  2‘  iV  3$  30 


(3-77) 


(3-78) 


p  =  Un2  +  P°2e 

21  cos  9' 


E22  =  ClVs2 


(3-79) 


(3-80) 


G21 


-C2E21"  CyOMl(A~  tanQ'^)’ClO  =  C3C4-G22  =  0  0 


(3-81) 


£21  =  (un2  +  Pu2)/c°s0'.^9,  =  C'.itanQ'  u9-  c^w^J 


21  5 


2  3  2e 


(3-82) 


H22=  "U2“  ^7  "2 


tort  0'  u;, 


(3-83) 


(3-84) 


T=c3cSw\ 


(3-85) 


C22  =  °>C21  =  C6(C7U2,  “  U'2)  +  C5  [C8U2  +  C3  U2  “  C9{W2e  +  <*N  “V 
Equations  (3-14)  through  (3-19)  are  already  in  a  form  similar  to  Equations  (3-70)  and  (3-78). 


(3-86) 


Finite-difference  schemes  will  be  developed  after  the  fin  boundary  scheme  is  developed. 
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3.7  BLUNT  BODY  SOLUTION 


As  earlier  indicated  in  the  geometry  description,  a  blunt  body  is  assumed  to  have  an  initial  circular 
cross  section;  however,  the  meridian  curve  need  not  be  a  circular  arc  segment.  At  the  nose  of  the  body,  the 
slope  angle  is  90  deg  and  is  larger  than  the  Mach  cone  angle.  Sharp  nose  bodies  with  initial  angles  greater 
than  the  Mach  angle  are  not  considered. 

The  basic  idea  is  to  start  the  surface  pressure  distribution  as  a  modified  Newtonian  one  and  continue 
the  computation  by  some  approximate  model.  This  approach  was  first  utilized  in  Reference  16  for  relatively 
large  supersonic  Mach  numbers.  The  second-order  shock  expansion  method17  was  utilized  to  continue  the 
computation. 

Two  different  methods  of  matching  a  Van  Dyke  second-order  axisymmetric  solution  with  a  Newtonian 
pressure  distribution  are  described  in  References  3, 18,  and  19.  Both  of  these  methods  use  surface  boundary 
conditions,  axial  source,  sinks,  and  doublets. 

The  first  thing  to  establish  is  the  effective  origin  of  disturbances  and  a  match  point  between  a  modified 
Newtonian  pressure  distribution  and  a  continuing  potential  solution.  In  Figure  5,  a  meridian  plane  and  the 
choice  of  match  point  and  origin  is  shown.  A  match  point  at  x  =  xm  is  chosen  so  that 


where  as 0.95  is  an  input.  xm  is 

A  tangent  from  m  is  drawn 
x  =  x0  S0.0. 


/dr  \ 

tan  8  =  —  =  a/(5,  a  <0.95 

m  \  dx 


(3-87) 


determined  by  Newton-Raphson  iteration. 

to  OT.  A  Mach  line  is  drawn  from  m  to  OU.  The  effective  origin  is  at  O', 


x  -  x  —  r,(x  ) 
o  m  b  m 


1  —  C 


dr. 

-  (x  ) 

dx  m 


+  cp 


(3-88) 


where  c  is  an  input.  Various  combinations  of  a  and  c  have  been  tried,  a  =  0.8,  c  =  0  seems  to  be  about 
optimum. 

Various  approaches  were  tried  for  a  model  blunted  cone  problem.  None  of  them  were  fully  successful. 
The  modified  Newtonian  pressure  coefficient  distribution  is  given  by 


C 


p 


P-P* 

y/2  Ml 


=  C  sin2 8 

po  e 


(3-89) 


sin  8  =  q  •  n. 
e  in 

=  cosa~sinS  -  sina~cos(Q-v)eos§ 

l  m  i  m 


(3-90) 


njn  is  the  inward  normal  to  the  body.  Cpo  is  the  stagnation  pressure  coefficient.  aT  is  the  total  angle  of  attack, 
and  v  is  the  roll  angle  of  the  crossflow  frec-strcam  velocity  component. 
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The  first-  and  second-order  potential  functions  are  given  by 

=  f0  ( x ,  r)  +  fx  (x,  r)  cos  (0  -  v  )  sin  aT 

<J>2  =  Xo  ('X'  ^  +  r) TOS  (Q  ~  v  )  sin  aT 
This  is  the  "hybrid”  approximation  of  References  5  and  19. 

The  potential  Bernoulli  pressure  coefficient  is  given  by 


(3-91) 

(3-92) 


C  = 

p 


1+  — -  M  2(1  -Q2) 
2  ® 


Y/(  y  —  1 ) 


-1 


Y/2  M 


(3-93) 


Equation  (3-92)  may  be  substituted  into  Equation  (3-93)  and  Equation  (3-90)  into  Equation  (3-89)  and 
matched  to  0  (sin3aT)  as  in  Reference  19.  Matching  coefficients  yields  Xox.  Xor.  fi,  fix.  and  firat  the  surface. 


It  is  then  assumed  that  the  ratio  of  first-  to  second-order  surface  velocities  at  zero  angle  of  attack  is  the 
same  as  the  conical  solution  ratio.  This  yields  an  approximation  for  fox  and  for  at  the  surface. 


The  remainder  of  the  initial  plane  of  data  is  obtained  by  utilizing  the  homogeneous  conical  potential 
functions  for  foxi  Con  Clx,  fir.  and  f^ 


In  summary,  an  approximate  initial  plane  of  velocity  data  is  obtained  that  matches  the  modified 
Newtonian  pressure  distribution.  The  solution  is  then  marched  downstream  by  an  explicit  routine.  Implicit 
routines  lead  to  strongly  oscillating  flow  prediction. 

The  blunt  body  option  should  only  be  used  for  significantly  blunted  bodies  that  have  large  wave  drags, 
since  the  CFL  step  size  (to  be  treated  later)  is  very  small  for  almost  pointed  bodies. 


3.8  FIN  EDGE  TREATMENT  AND  SURFACE  DISCONTINUITIES 

The  first  objective  is  to  obtain  jump  conditions  for  supersonic  leading  and  trailing  edges.  For  a  fin 
leading  edge,  the  velocity  ahead  of  the  edge  is  known  and  the  velocity  on  both  edge  sides  is  to  be  determined. 
For  a  variable  sweep  leading  edge,  the  envelope  of  disturbances  or  wavefront  from  the  edge  is  rotated  with 
respect  to  the  local  0  plane  of  the  lifting  surface. 

Conservation  relations  for  the  first-order  wave  equation  may  again  be  derived  by  the  method  of  "weak 
solutions”. 


(AtPj)2  +  (Atfj)2  -  (32(AUj)2  =  0 

Here,  the  conservation  jump  relations  for  the  three  marching  equations  have  been  combined. 

At  k  =  l,  the  remaining  relationships  are  determined  by  satisfying  both  the  body  and  fin  boundary 
conditions  where  the  fin  meets  the  body  or  the  beginning  of  the  corner.  At  k  *  1 ,  the  conservation  of  the 
tangential  velocity  component  and  the  fin  boundary  condition  provides  the  needed  two  relations  for  closure. 

At  the  fin  surface  (from  Figure  2), 
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tp  =  9  -  8^.  -  o=  0 


(3-95) 


Equation  (3-95)  may  be  used  to  determine  the  fin  surface  normal.  The  fin  full  boundary  condition  may 
be  written  as 


W  =  ro  U .  +  ro  V, 

1  i  1  r  l 


(3-96) 


ro  =  rd/c tic 

X 


tan  ~^(Ur) 


1  +  (t/rf 


~  t  =  tan 8.  ±  t 
x  f  ox 


tr-tir 

ro  —  — -  ~t  —  tlr 

r  .  .  v2  r 


1  +(t/r¥ 

Here  8f  is  the  fin  control  deflection,  tox  is  the  fin  slope;  ( + )  for  0  >  0f  and  ( — )  for  9  <  9f  sides  of  the  fin. 
The  fin  boundary  condition  is  further  approximated  by  neglecting  the  spanwise  contribution. 


W,.+Au),  =  t 

16  1  x 


(l-Fl^  +  FO^  +  AUj) 


(3-97) 


(3-98) 


In  almost  all  applications,  F  =  0  or  the  approximate  fin  boundary  condition  yields  a  result  superior  to  that  for 
F  =  1. 


For  k*  1,  the  tangential  velocity  may  be  written  as 


V  =  sin  \ 


U  cos  8  ^  +  W  sin  8^. 


(3-99) 


+  V  cos  \ 

where  \  is  the  local  sweepback  angle.  Applying  conservation  of  the  tangential  velocity  component  yields 

(3-100) 

At’j  =  -  tanACAUjCosS^.  +  A w^sindj.) 

In  all  cases,  a  quadratic  equation  must  be  solved  for  the  jump  velocities.  Choice  of  sign  is  determined  by 
the  turning  angle  being  compressive  or  expansive. 

The  case  for  k  *  1 ,  constant  Uib  =  1 ,  Wjb  =  -  sin  a,  8f  =  0,  tox  =  0  may  be  readily  solved  and  is 
illustrative. 


Aui=*-r 


stria 


W-tan2\)m 


(3-101) 


•  /  %»  r*  * ~,h  •  r  *  -  m  f  ' 
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This  is  the  well-known  result  for  the  constant  supersonic  region  for  a  thin  wing  due  to  Ackeret.  The  solution 
indicates  that  as  X  increases,  the  jumps  in  velocity  gets  very  large  and  that  small  disturbance  theory  is 
violated  locally.  The  subsonic  leading  edge  case,  for  which  conservation  relations  do  not  exist,  occurs  when 
tan  X  >  p.  Such  singularities  do  not  arise  for  the  full  nonlinear  potential  equations. 

A  maximum  sine  of  the  sweepback  angle,  sin  Xm,  is  defined  by  requiring  that  the  Mach  wave  angle  be 
swept  back  at  least  X0  from  the  leading  edge. 


sinX  =  (  —  sin  X  +  pcos\)/M 
m  o  o  °° 

where  sin  X0  is  a  code  input.  For  X  >  Xm,  a  subsonic  leading  edge  heuristic  treatment  is  used. 

In  order  to  evaluate  the  nonconservative  gradient  jump  contributions,  it  is  assumed  that  the  local  wave 
front  orientation  for  the  second-order  is  the  same  as  for  the  first-order  case. 

The  wave  front  angle,  c,  is  defined  clockwise  from  the  fin  plane.  On  a  wave  front,  the  tangential 
component  of  velocity  is  conserved  in  the  local  r,  0  plane,  which  yields 


tone  = 


A  w 

l 


(3-103) 


For  the  case  of  Equation  (3-101), 


_  fanX  _  tanX 

tone  =  +  — - - — — ,  sine  =  H - — 

.2^V2’  p 


(3-104) 


(P*-  ton  X) 

This  is  also  a  well-known  result  from  conical  flow  theory  of  a  fiat  plate  swept  wing. 

A  second-order  application  of  the  method  of  "weak  solutions"  and  Equation  (3-103)  yields 


(Aio2)2  +  (Ao^2  -  p2  (Attg)2  =  p2  A u'  Au2  +  Ai/  Au2  +  Aw'  A u>2 


(3-105) 


M 

Au'  =  —■  AA 
p2 


(3-106) 


Ay' 


2  — Y 


AT 


— -  =  A  q  —  -  — —  A  v 

■'2  2  r2 


(3-107) 


M 


Aw' 
M 2 


=  a70- 


2-y 


M 


6  2  p2 


A  w 


(3-108) 
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Au  =  “  A(Oj)  +  A(WjL»t)  +  2/3  tone  A(u)j)  -  (ulfe  +  u)lh tone) 

A  w  =  -  A(wf)  +  A(v2  w  )  4-  2/3  co/eA(v?)  -  (w..  4-  u,,  cote)  A(y^) 

31  II  1  lO  lO  1 

Closure  equations  comparable  to  (3-98)  and  (3-100)  are 


(3-109) 


(3-110) 


W2fc  +  Au,2 


t 

X 


(1  -F)Ula  +  F(U2b+*u2) 


Ao 


2 


—  ton  X 


A  u  2  cos  8^.  +  &w2sin8j. 


(3-111) 


The  supersonic  trailing  edge  requires  different  treatment.  Here,  the  flow  goes  from  upper  and  lower 
surfaces  to  the  free  stream.  A  formal  full  analysis  would  require  an  iterative  solution  for  the  flow  turning 
angle.  In  many  cases,  no  solution  would  be  possible.  The  simplest  approach  is  to  assume  that  the 
downstream  velocity  vector  lies  in  the  plane  of  the  fin.  This  is  the  same  as  tx  =  0. 


At  the  trailing  edge  ’  there  is  a  velocity  and  pressure  jump  on  both  sides  of  the  arbitrary  wake  cut  (the  fin 
plane).  The  resultant  pressure  coefficients  and  v  component  of  velocity  are  averaged  across  the  wake.  At  the 
corner,  the  v  velocity  component  is  not  averaged  and  the  body  boundary  condition  is  used. 


Turning  angle  for  the  corner  at  k  =  1  is 


DTH  =  -sin-1 


-^26±(-to"S^+VV 


'[kj 


u+fV72 


Turning  angle  away  from  the  corner  for  k^  1  is 


(3-112) 


DTH-  -sin-1 


-^^26-^±(-^5ft/26+Wr26) 


(Gj“  Va+V 


1/2 


(3-113) 


Equation  (3-101)  illustrates  the  problem  for  subsonic  leading  edges  or  edges  that  are  slightly  supersonic. 
The  quadratic  solution  for  the  velocity  jump  is  imaginary  or  very  large.  For  a  delta  wing  in  a  free  stream,  the 
first-order  conical  solution  due  to  Busemann  is  well  known.  The  pressure  coefficient  has  a  square  root 
singularity  near  the  leading  edge.  The  subsonic  leading  edge  pressure  distribution  is  assumed  to  behave  like 


,k*l 


(3-114) 


Here,  is  the  value  of  I,  at  the  leading  edge  and  ^  is  the  value  of  £  at  a  point  near  the  leading  edge.  Inner 
edges  are  always  assumed  to  be  supersonic.  The  root  point,  k  =  1 ,  is  always  supersonic.  Ajj  is  evaluated  by 
setting 


*-l 
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Alternates  to  the  expression  (3-115)  were  tried,  but  proved  to  be  somewhat  less  successful.  One  of  these 
alternatives  was  the  mean  value  relationship. 


S(*-D 


K1 


C  d  = 
p  t 


\*-i> 


SKJ 


«J  -  S2)172 


di. 


Here,  KI  is  the  value  of  k  near  the  inner  edge. 

The  above  relations  hold  for  first-  or  second-order  and  either  side  of  the  fin. 


(3-116) 


Closure  conditions  for  the  subsonic  edge  were  given  by  the  requirement  that  the  velocity  vector 
component  parallel  to  the  edge  is  conserved  and  the  downstream  boundary  condition  is  satisfied.  The  former 
requirement  is  given  by  Equation  (3-99). 

The  boundary  conditions  are  to  first-  and  second-order: 


W.  =  t 

1  X 


(1  -F)u  +FU , 


W=t 

2  x 


(1  -F)U1  +  FU2j 


The  velocity  vector  squared  is  determined  from  the  value  of  Cpk 


(3-117) 


U2  +  V2  +  W2  =  1  - 


Cpk}K 


+  1 


Y-l 

Y 

-1 


(3-118) 


At  a  subsonic  side  edge,  it  is  assumed  that  the  value  of  the  velocities  just  leaving  the  edge  are  equal  to 
value  at  next  k  point  out;  i.e.,u it  =  (ui)  k+i,  etc. 


3.9  FIN  SURFACE  ADVANCEMENT  EQUATIONS 

Characteristic  compatibility  relations  for  the  fin  advancement  equations  are  sought  for  the  qr,  x  plane 
with  £  as  an  interior  variable  at  ip  =  0. 

Starting  equations  for  first-  and  second-order  combined  are 
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l 

u  —  u  a  +  —  a  v 
x  ip  j  p2  r  > 
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1 


rp2  W'V  P2 


1  1 

y/r+y^-H-u,^  — 


aA  l 

'  ,  %  l 

^9  ^9  11  qN 

o  —  +  - 
*  dip  r 

d/diyirq)^  —  a/5  tp  (r  qj  +  — 

N* +i^ 

aA  8A 

4-  £  — — 

ax  *  a$ 


-  ^d/d^ruj 


ar  a(rUj)  j 

r  aiy  r 


dwl  3  u>j 

L  — -  4-  - 

^  a0  avp 


-u  A 

1  x 


(3-119) 


au 

y  — o  y  -fo  —  =  -  y.  E,  +  ur£ 
*  *  ^  r  aiy  i  *  t  > 


(3-120) 


ty  — 

X 


1  1 
aw  —  -  a  =  —  iyt  £,  H —  u.L 
i  ip  r  4<  lx  r  1  v 


(3-121) 


Solution  for  the  eigenvalues  and  vectors  are  similar  to  that  for  the  body 

[l+(ro  )¥2 


A  =  -o  + 

1  X 


Pr 


•X2=  -°x 


(3-122) 


Choice  of  sign  is  dictated  by  the  incoming  characteristic  on  either  side  of  a  fin.  The  first  sign  is  for  the 
0  >0fside. 


Left  eigenvectors  are 


L]/L\  =  ±PU+(nir)2]w 


=  l,ij  =  ror 


(3-123) 

(3-124) 


The  surface  fin  boundary  advancement  equations  are  evaluated  at  ip  =  0,  where  d/dqi  ->a/30  for  the  thin- 
fin  approximation.  In  the  second  advancement  equation,  0  may  be  eliminated  by  using  the  unused  third 
irrotationality  equation  as  was  done  for  the  body. 


Terms  involving  ror  are  dropped.  However, 


dJdx  {rv  )  -  —t  I  r 
r  x 


(3-125) 


is  retained.  The  first-order  advancement  equations  may  be  written  as 


fn  +fi!  +  fi=c 

ax  di,  <30  11 


(3-126) 
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1 

t 

I 

dEl2  ^12 

—  + 1  - - - 1  =  c 

toe  *  0$  a$  12 

(3-127) 

i 

£&  u  +  w, 

E„=— i— > 

(3-128) 

' 

C4 

C4  has  been  defined  earlier  in  Section  3.6. 

E12  =  Ul/c4 

(3-129) 

„  w  V, 

Gn~~  et/r  («!±  j)T  j+Es(±pUl  +  Wl) 

(3-130) 

Et  and  Ex  have  been  defined  in  Section  3.3. 

G12  =  wi  ~  vxEt/R 

(3-131) 

Hl  =  ~dl^tx(iPui  +  u;i)+  Uj  ±  IOj/PUj  =  (1  /P-Rb)IR 

(3-132) 

! 

Ch  =  <W±^',>-<M±  J 

(3-133) 

C12  =  d2txVl 

(3-134) 

The  second-order  advancement  equations  are  of  the  form 

+  ,as  w 

toe  %  00  21  +  +  a0) 

(3-135) 

o8 

II 

si'* 

(3-136) 

£21  ~  (±P  u2e  +  WJ,C4 

(3-137) 

E22  =  U2  /c4 

(3-138) 

ET  V2 

G21  =  “  T  (U2  *  *  J  +  *  (±P  +  ifj) 

(3-139) 
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^22  W2  ETV2*R 


H=  -d. 


tx{±VU2e  +  W2)+  U2±W2e,V 


C2X  =  <*1 


,  1 

,1 

±  ~  (o„  +  q„v.)  +  t  v~ 
n  2e  ^  Jv  1  i  2 

d2 

±Viw2'  +  (lNWl)+U2 

C22= ' d2  V V S =  ±  J  (W  1  +  ET  Wx'R)’  T=  ±  dX  »l'P 


(3-140) 

(3-141) 

(3-142) 

(3-143) 


3.10  FIN-BODY  CORNER  ADVANCEMENT  EQUATIONS 

The  fin-body  junction  is  a  crossflow  geometry  discontinuity.  In  Reference  10,  the  corner  is  treated  as  a 
single  valued  point  or  streamline. 


A  full  characteristic  analysis  involves  the  matrix  equation  form 


P  +  A  P.+  BP  =R 

x  l,  V 


(3-144) 


The  body  and  fin  analysis  have  already  considered  characteristics  in  the  S  =  0  and  ip  =  0  planes.  In  the  current 
case,  S=0,  ip  =  0  at  the  corner.  In  Reference  10,  no  unique  solution  was  found,  if  one  considers  the  corner  as  a 
streamline  and  hence  isentropic  for  the  Euler  equations.  In  Reference  10,  a  heuristic  analysis  led  to  two 
alternate  differential  equations  for  advancing  p.  p  combined  with  isentropy,  total  enthalpy  and  the  corner 
vector  is  enough  to  advance  all  variables. 


Here,  a  nonunique  set  of  relations  is  established  by  advancing  Q2/2  and  using  both  boundary  conditions 
for  a  streamline  along  a  corner 


1  dQ2 

2  dx 


l 

- -a/ax 

(y-VmJ 


(3-145) 


drb 

V=  —  U  +  tanQ'W 
ax 


(3-146) 


W=tU 

X 


(3-147) 


Equation  (3-145)  may  be  combined  with  the  three  vector  velocity  advancement  equations.  A  combined 
first-  and  second-order  nonconservative  advancement  equation  for  (3-145)  may  be  written  as 


a£  dE 

—  -M  — 

dx  x  a$ 


t  X  dll 
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r  r  *  r  °  3^  r 
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1  1  au>i  1  flu>. 
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(3-148) 


All  the  transformation  derivatives  are  to  be  evaluated  at  £  =  0, 0  =  0r.  In  addition, 


U  =  U  + 


U  2  +  V2  +  W2 

e 

~~  2 


For  First  order,  A  =  0,  ve  =  vj ,  and  we  =  w* .  For  second  order,  ve  =  V2e  and  we  =  W2e. 


(3-149) 


3.11  EXPLICIT  FINITE-DIFFERENCE  SCHEME 

The  MacCormack  scheme20  is  applied  to  body  surface  points,  Field  points.  Fin  surface  points,  and  Fin-body 
corner  points.  For  all  these  points,  the  E  vectors  are  to  be  advanced  by  a  predictor-corrector  Finite-difference 
scheme.  For  the  field,  there  are  three  vector  equations;  for  the  fin  and  body,  there  are  two  vector  equations; 
and  for  the  corner,  there  is  a  single  equation. 


E*  is  computed  utilizing  (3E/3x)i  which  is  evaluated  by  using  forward  or  backward  one-side  differences  for  £ 
and  9  derivatives.  E*  is  utilized  to  decode  for  the  new  velocity  components  at  the  predicted  conditions.  The 
(3E/3x)*  vector  is  utilized  to  obtain  the  final  trapezoidal  rule  value  of  E;  + 1.  The  reverse  backward  or  forward 
one-sided  differences  for  E,  and  9  derivatives  is  used  to  evaluate  (3E/3x)*  at  the  predicted  state.  A  decoding  of 
the  E;  + 1  vectors  yields  the  advanced  velocities.  For  field  points,  a  code  option  uses  forward  E,  and  9  differences 
for  the  predictor  and  backward  £  and  9  differences  for  the  corrector  or  alternate  order  from  step  to  step.  At  the 
surfaces,  the  boundary  conditions  are  also  used  in  the  decoding. 

At  a  solid  body  surface,  only  forward  differences  in  $  are  possible.  If  k  =  2  lies  on  a  fin  point,  the  E, 
derivative  is  set  to  zero. 

For  fins,  the  0  differences  must  be  always  forward  on  the  9  >  0f  side  and  backward  on  the  0  <  9f  side. 
Points  on  a  fin  next  to  an  edge  need  special  consideration.  For  supersonic  edges,  £  differences  for  points  like  A 
or  B  in  Figure  2  are  always  backwards.  For  similar  points  on  an  inner  edge,  5  differences  are  forward.  For 
subsonic  outer  edges,  a  forward  difference  to  a  point  in  the  flow  such  as  C  is  taken.  For  the  point  C  in  the  flow, 
the  £  difference  for  a  supersonic  edge  is  always  forward.  For  a  subsonic  edge,  a  backward  difference  to  the 
average  of  A  and  B  values  is  utilized. 

For  explicit  hyperbolic  schemes,  a  CFL  condition  for  stability  is  required.  A  formal  analysis  was  not 
made.  The  step  size  is  chosen  so  that  the  upstream  Mach  cone  from  a  point  at  x  4-  Ax  lies  within  the  nearest 
grid  points  of  the  plane  at  x. 

The  radius  of  the  Mach  cone  at  x,  0  +  A0,  and  E,  +  A$  to  first  order  is 

i 
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drb 

—  (i-OAe  +  *a/p  -Rb)^ 


2 


For  a  constant  £  value,  Equation  (3-151)  is  a  minimum  at  the  surface. 


(3-151) 


For  a  constant  0  ray, 


(3-152) 


(3-153) 


Equation  (3-153)  is  a  minimum  at  £  =  0  as  well  in  most  applications.  A  search  as  8  varies  for  Equations 
(3-152)  and  (3-153)  is  made  for  (Ax)  min.  The  step  size  chosen  is  0.9  of  this  value. 


Discontinuities  create  large  differences  in  advancement  point  velocity  values.  These  differences  must 
somehow  be  dissipated  for  the  First  few  steps  after  the  discontinuity. 


For  a  body  planar  discontinuity  (e.g.,  a  slope  or  inlet  discontinuity),  the  (j  derivatives  at  (;= 0  are  set 
equal  to  zero  for  the  first  step  and  increased  linearly  to  the  full  value  after  i^  steps,  which  is  a  code  input.  A 
curvature  discontinuity  does  not  have  to  be  damped.  However,  much  better  results  are  obtained  by  marching 
up  to  the  curvature  discontinuity  rather  than  across  it.  An  ij  of  3  or  4  steps  is  usually  adequate  for  damping 
body  discontinuity  disturbances. 


At  a  Fin  leading  edge,  the  0  discontinuity  disturbance  can  be  very  large  and,  for  swept  back  Fins,  it  can 
continually  be  created  until  all  (j  traces  cross  the  leading  edge.  Both  ^  and  0  derivatives  can  be  severely 
affected,  in  particular  for  a  subsonic  leading  edge.  The  root  velocities  are  also  affected. 

Only  damping  of  the  ©derivatives  just  past  the  leading  edge  is  used.  The  number  of  steps,  ijf,,  for 
damping  is  computed  from  the  turning  angle  of  Equations  (3-112)  or  (3-113). 


V  =  INT 


\9.55  DTH\  +  0.5 


(3-154) 


INT  is  the  integer  conversion  of  the  floating  point  value.  For  a  subsonic  edge,  the  number  of  steps  is  doubled. 
At  the  root,  ijf  =  6  is  chosen  for  a  subsonic  edge. 


Damping  of  8  derivatives  for  the  Fin  is  similar  to  £  derivative  damping  for  a  planar  discontinuity.  The 
derivatives  are  linearly  increased  from  zero  to  the  Final  value  after  ijf  steps.  At  the  root,  no  damping  is  used 
fc”  the  first  ijf  steps.  The  velocity  vectors  are  advanced  using  the  corner  jump  relations  twice:  for  the  first 
jump,  the  body  shape  is  advanced;  and  for  the  second  jump,  the  fin  shape  is  advanced. 


Due  to  the  strong  velocity  oscillations  on  a  fin  with  subsonic  leading  edges  and  for  subsonic  ridge  or 
surface  discontinuity  lines,  no  attempt  was  made  to  compute  jumps  for  fin  ridge  discontinuities.  For  a  fin 
with  all  supersonic  edges,  the  jumps  could  be  computed. 


It  is  assumed  that  for  a  good  design,  the  body  drag  contributes  most  of  the  wave  drag  and,  hence,  the  fin 
wave  drag  may  be  approximated.  The  normal  force  contribution  of  a  fin  is  given  fairly  accurately  by  the  plan- 
form.  The  wave  drag  of  the  fins  is  given  by  the  thickness  distribution  and  deflection  angle  of  the  fins.  The 
characteristic  compatability  boundary  condition  will  capture  small  jumps  on  the  fin  surface.  An  equivalent 
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thickness  distribution  would  bo  given  by  multiplying  the  actual  tox  by  a  thickness  reduction  factor,  trf,  and 
then  accounting  for  the  actual  fin  thickness  when  the  axial  force  contributions  are  computed.  This  will  be 

elaborated  on  later. 

In  order  to  dampen,  somewhat,  the  very  large  amplitude  oscillations  in  the  surface  pressure  for  lifting 
surfaces  with  subsonic  or  nearly  subsonic  edges,  a  smoothing  operation  was  introduced  for  the  E  vectors 
before  decoding  E(  ¥\. 

A  weighting  function  is  chosen  as  Qs  =  Vi2  +  Vi2  +  Wi2,  which  is  a  positive  definite  quantity.  Field, 
body  surface,  fin  surface,  and  corner  vectors,  Ec,  are  smoothed  as  follows: 

Field 


Both 


E  =  (1-C.  -C  -C.  .  -  C.  u.)£c.. 

jM  j- 1,*  7  +  1,#  7,*-l  J,  k+ 1  j,k 

+  C.  . .  Ef.  .  .  +C...  ,  E°  .  +C.  .  .  E°.  .  . 

j-l,k  j-l,k  7  +  1,*  7  +  1,*  7,*-l  7,*-l 

+  C  .  .  ,  Ef.  ,  , 

7,  *+l  7,  *  + 1 


(3-155a) 


Fin 


Corner 


E.=(l-C.  .  .  -  C.xl  .)EC  ,  +  C  Ef.  +C  ,  E0  , 

7.1  7-1.1  7  +  1,1  7,1  7-1.1  7-1,1  7+1,1  7+1,1 


(3-1 55b) 


E  =(1— C.,  ,-C.  li))£c  .  +  C.  .  ,EC,  ,+C..i1Ec.i1 

J,k  7,*-l  7,  k+ 1  j,k  7,  *—  1  7,  *  —  1  7,fe+l  7,*+! 


(3-155C) 


Here, 


ej  ,k  =  <l-cj±x.rcjJEf>j.i  +  c7±i.i^±u  +  S.2^.2 


(3-155d) 


C.+1  .  =  C 

7±1,*  s 


«Ajk  +  «A-±i.* 


(3-156a) 


Cj.k±l 


Wj.k-<Vj.k±i 


^s)j,k+^s)J,k±i 


(3-156b) 


C8  is  a  smoothing  constant  level,  typically  0.2. 

Note  that  on  the  0>0f  side  of  the  Nth  fin,  the  data  are  stored  at  j  =  JM  +  N. 
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The  damping  smoother  does  somewhat  limit  the  amplitude  of  oscillations.  The  oscillations  die  down  far 
from  the  leading  edge. 

The  explicit  marching  solution  will  fail  when  the  value 

1  +  ^-  M2(l-Q2) 

2  00 

becomes  negative  in  Equation  (3-93).  To  avoid  a  complete  halt,  the  Bernoulli  quantity  is  set  to  zero  at  a  solid 
surface  and  the  pressure  becomes  the  vacuum  value.  Readjustment  of  the  total  velocity  components  requires 
two  additional  relations.  One  of  these  relations  is  given  by  the  fin  or  body  boundary  conditions:  at  the  corner, 
both  boundary  conditions  are  used;  on  the  body,  the  W  component  remains  unchanged;  and  on  the  fin,  the  V 
component  remains  unchanged.  Sometimes  useful  results  can  be  obtained  by  this  delay  of  solution 
breakdown.  A  complete  breakdown  in  the  computation  occurs  when  the  vacuum  values  of  Cp  are  propagated 
along  the  solid  surfaces. 


4.0  HIGH  MACH  NUMBER  INVISCID  SOLUTION 

At  high  Mach  numbers,  the  origin  Mach  cone  crosses  pointed  body  surfaces  and  other  compressibility 
effects  become  dominant.  Here,  an  adaption  of  the  local  pressure  solution  methods  of  Reference  9  is  used.  The 
local  pressures  depend  primarily  on  the  orientation  of  the  local  normal  with  respect  to  the  free-stream 
velocity  vector.  These  methods  are  extended  and  applied  in  References  21  and  22.  The  methods  are 
summarized  in  Table  1.  The  methods  use  Newtonian,  tangent  cone,  tangent  wedge,  and  Prandtl-Meyer 
pressures  and  empirical  blends  of  these. 

As  indicated  in  Section  1.0,  Euler  code  solutions  such  as  described  in  Reference  10  give  good  solutions  for 
complex  configurations.  The  local  solution  methods  run  much  faster  than  the  marching  Euler  solution  and 
provide  rapid  trade-off  computations;  however,  they  sacrifice  considerable  accuracy. 


5.0  INVISCID  LOADING  COEFFICIENTS 


Axial  Force  Coefficient 


AR 


■2" 

o  r>>  ax 


SF 

<Ve  +  I 


n  =  1 


(C  t  -c 

ppn  xpn  pmn 


t 

xmn 


)  dr 


(5-1) 


Here:  n  refers  to  the  fin  number;  Cpp  is  the  pressure  coefficient  on  the  0  >  8f  side;  Cpm  is  the  pressure 
coefficient  on  the  0<9f  side;  txp  and  txm  are  corresponding  body  slopes;  rj  and  r0  are  inner  and  outer  fin  edge 
radii;  and  b  refers  to  the  body. 


Due  to  large  control  deflections,  fin  discontinuities,  or  subsonic  leading  edges,  some  finned  body 
configurations  cannot  be  computed  at  the  given  geometry.  The  computation  is  made  at  reduced  values  of  toxn 
=  txon.  The  adjustment  to  the  fin  contribution  is  given  as 


C  t  -C  t  =  (C  -  C  )tan8,  +t  It2,  (C  +C  ) 

ppn  xpn  pmn  xmn  ppn  pmn  fn  xon  rfn  ppn  pmn 


(5-2) 
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Prime,  here,  refers  to  the  reduced  slope  geometry  values.  One  uncertainty  is  whether  or  not  the  second  term 
m  Equation  (5-2)  will  be  lost  in  the  noise  of  the  oscillating  pressure  contribution  to  lift.  Here,  trfn  =  txon/toxn. 

Normal  Force  Coefficient 


Pitching  Moment  Coefficient 


dC 
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*  dx 


2n  3rfe 

(r.cosQ  +  sind)C  .dQ 

a  5  ae  p*> 


NF 
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n  =  1 
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(C  —  C  )sm8.  dr 

ppn  pmn  frx 


(5-3) 


dC  dCKt 
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(5-4) 


+  I 

n=  1 


on 


r. 

in 


rcosQ.  (C  t  —  C  t  )dr 

fn  ppn  xpn  pmn  xmn 


It  is  not  certain  whether  the  adjustments  of  Equation  (5-2)  are  really  needed  for  the  finned  contribution  of 
moment  due  to  drag. 


Large  combinations  of  a  and  5  do  not  permit  inviscid  computation  of  pitching  moment  and  normal  force 
in  the  pitch  plane.  Computation  at  reduced  values  will  permit  extracting  the  linear  dependences.  First,  the 
computation  is  made  at  a'  with  no  control  deflection.  Then,  the  computation  is  made  at  a'  with  a  control 
deflection  5'.  This  is  readily  done,  since  the  computational  code  saves  computational  values  up  to  a  section 
end  and  can  be  restarted  at  the  section  end  and  run  to  the  end  of  the  body.  For  more  than  one  set  of  control 
deflections,  additional  computations  may  have  to  be  made. 


The  adjusted  inviscid  normal  force  coefficient  is  then  given  by 

CKJ  =  C.r  sin  a/sin  a'  +  AC. ,  ton  8/ tan  5' 

/V  |V  IV 


(5-5) 


AC'n  is  given  by 


ACn  =  CN(a',8'*o)-  CN(a\ 8  =  0) 


(5-6) 


The  pitching  moment  coefficient  is  treated  in  a  similar  manner.  The  remaining  coefficients  are  assumed  to  be 
computed  at  small  angles  of  attack  and  side  slip.  For  configurations  without  quarter-plane  symmetry,  a  and 
a'  must  be  relative  to  the  Cm  =  0  angle  of  attack.  Note  that  Cn  and  Cm  do  not  go  to  zero  simultaneously  for 
half-plane  symmetry. 
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Side  Force  Coefficient 
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Rolling  Moment 


LRARCt 


SF 

+  I 

n  =  1 


f2" 

Jo 


C  Kd0 

pb 


r. 

in 


r(C  -C  )  dr 
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(5-9) 


The  0  integrals  need  not  be  taken  to  2n  in  all  cases,  except  for  LM  =  6.  For  no  fins,  the  0  integration  is 
continuous.  For  a  finned  body  section,  the  0  integrals  are  piecewise  continuous.  0  integration  is 
approximated  by  a  three-  or  four-point  Simpson’s  rule.  This  implies  a  minimum  of  three  points  between  fins. 
The  0  variable  step  integration  is  converted  to  <}>  equal  step  integration  by  the  80/8<J>  Ac|>  values. 


The  r  integrations  are  carried  out  with  the  aid  of  the  $  and  £  transformations.  A  Simpson’s  three-point  or 
trapezoidal  rule  integration  is  used  for  the  points  between  k=  KI,  KO.  KI  and  KO  do  not  lie  on  the  edges.  The 
end  interval  integrations  utilize  an  extrapolated  trapezoidal  rule. 


After  the  body  section  computation  is  completed  ,  the  coefficients  are  integrated  from  the  origin  to  a 
value  of  x(i.e.,  a  load  coefficient  buildup).  For  a  section  with  no  fins,  a  piecewise  leap-frog  three-point  or 
quadratic  approximation  for  the  coefficient  gradients  is  utilized.  For  a  section  with  fins,  a  piecewise  linear 
approximation  is  utilized  for  the  coefficient  gradients.  Note  that  those  integrations  are  not  the  same  as 
trapezoidal  and  Simpson’s  rules. 


6.0  SKIN  FRICTION,  BASE  PRESSURE,  AND  HIGH 
ANGLE-OF-ATTACK  PREDICTION 


6.1  SKIN  FRICTION 

The  approach  here  is  that  described  in  Reference  23.  The  body  skin  friction  is  based  on  the  body  surface 
area  and  body  length  Reynolds  number.  The  fin  skin  friction  is  based  on  twice  the  fin  planform  area  and  the 
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mean  aerodynamic  chord  Reynolds  number.  The  mean  skin  friction  coefficient  is  given  by  the  Van  Driest 
model  of  Reference  24. 

The  skin  friction  coefficient  for  a  body  or  fin  element  is  given  by 


C - - 

AF  Ai 


Rnc  (  1.328 


CFL+  R  °FC) 

NL  v  n  NC 


f°rRNL>RNC 


(6- la) 


Sw  1,328 
AR  VRNL 


f°rRNL<RNC 


(6-lb) 


Here  Sw  is  the  element  wetted  area,  Ar  is  the  reference  area,  Rnc  is  the  critical  Reynolds  number,  and  Rnl  is 
the  element  Reynolds  number.  Rnc  f°r  the  body  is  chosen  as  1.0  xlO6  and  0.5  x  106  for  a  fin.  The  Reynolds 
number  is  given  by 


(6-2) 


R  mfm  is  the  Reynolds  number,  per  Mach  number,  per  foot,  per  million,  which  is  a  function  of  altitude  for  a 
given  atmosphere.  Lc  is  a  conversion  factor  for  converting  the  element  length,  L,  to  feet.  Crl  and  Cpc  are 
mean  skin  friction  coefficients  at  the  Reynolds  numbers  Rnl  and  Rnc.  respectively. 


The  mean  skin  friction  coefficient  is  given  by 

dllVCF-bgCF  =  d2 

where 


(6-3) 
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Cf  in  Equation  (6-3)  is  solved  by  iteration  at  Rn  =  Rnl  and  Rn  =  Rnc- 

L  for  the  body  is  the  body  length.  L  =  Lf,  the  mean  chord,  for  the  fin  for  a  nontrapezoidal  planform  is 
estimated  for  an  equivalent  trapezoidal  planform. 
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From  Figure  3, 


rXMA 

isw  =  Af=  *SMA(C  +  c,  )  =  (SO—SDdx 

rm  tm  iXMl 


(6-4) 
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j  XMA 

-  SMA2(C  +2C.  )  =  i  (SO2  -  SI2)dx 

6  rm  lm  XMI 


(6-5) 


Crm  and  Ctm  are  the  equivalent  trapezoidal  planform  root  and  tip  chords,  respectively.  Use  Equations 
(6-4)  and  (6-5)  to  solve  for  Crm  and  Ctm  and,  then, 


L  =  2  <■  -  <6'6> 
f  3  (C2  -  C2  ) 

rm  lm 

This  approach  avoids  having  to  compute  the  chord  at  the  spanwise  location  of  center  of  planform  area. 

6.2  BODY  BASE  PRESSURE  AXIAL  FORCE 

A  solid  noncircular  base  is  considered.  No  extensive  set  of  data  bases  seems  to  be  available. 
Approximate  wake  boundary  layer  interaction  models  or  Navier-Stokes  model  solutions  are  not  compatible 
within  the  current  context  of  rapid  computation. 

A  rational  approximation  is  developed  that  utilizes  a  minimum  of  available  data: 


C  (0)-  C 

pr  po 


r2/r2  (0) 


AT  (0) 


K+  C,  (AT)  (1-/0 

pZ 


(6-7) 

(6-8) 


Cprr  is  the  base  pressure  coefficient  assuming  a  local  free-stream  Mach  number,  M  -  In  Reference  26, 
Chapman  correlates  base  pressure  data  using  Pb/p=°  =  Pb/p  p’/p®  for  axisymmetric  bodies.  p'/p«,  is  the 
pressure  distribution  at  the  end  of  a  body  extended  by  a  pseudo  cylinder.  The  cylinder  has  a  diameter  and 
length  equal  to  the  base  diameter.  For  a  noncircular  cylinder,  the  cylinder  length  may  be  assumed  to  be 
equal  to  the  equivalent  circular  diameter  based  on  the  base  area,  (4  At>/n)L  This  approach  is  possible  for  the 
low  Mach  number  regime.  The  base  pressure  distribution  then  becomes 


C  =  C  (AT )2  (1/M2  +  0  .7  C  )  -  C 

pr  prr  °°  p  p 

For  higher  Mach  numbers,  this  approach  does  not  seem  feasible  and  Cpr  is  assumed  to  be  Cpr  =  Cprr. 


(6-8a) 
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Here  M'  (8)  is  the  local  Mach  number  at  the  end  of  the  body,  Cp2  (M1)  and  CP3  (M')  are  the  mean  base 
pressures  from  References  25  and  26,  respectively.  CP2  is  the  base  pressure  for  a  wide  and  long  plate  and  CP3 
is  the  base  pressure  at  the  end  of  a  long  cylinder,  k  is  the  local  curvature 


k(0)  =  - 


(6-10) 


Kma  is  the  maximum  curvature.  For  a  circular  base,  K  =  1  and  Cprr  =  Cp3.  However,  for  a  circular 
base,  there  is  still  an  angle-of-attack  dependence.  The  Cp  variations  in  Equations  (6-7),  (6-8),  and  (6-8a)  are 
clearly  arbitrary  and  need  some  correlation  with  data. 


The  base  axial  force  coefficient  is  given  by 
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(6-11) 


At  low  Mach  numbers,  the  Cpr  computation  given  by  Equation  (6-8a)  is  sensitive  to  the  M'  computation 
for  bodies  with  lifting  surfaces.  In  this  case,  the  value  of  Cab  given  by  Equation  (6-8)  should  also  be  computed 
with  no  body  extension.  A  minimum  of  the  two  computations  is  the  best  choice. 


6.3  HIGH  ANGLE-OF-ATTACK  PREDICTION 

The  methods  of  Reference  2  predict  nonlinear  charcteristics  for  the  pitch  plane  as  well  as  rolled 
orientation.  These  methods  combine  nose  vortex  tracking,  surface  singularity  solution,  leading  and  side  edge 
suction,  and  fin  shed  vortex  tracking. 

An  alternate  approach  is  to  adapt  the  crossflow  lift  extension  of  Reference  27. 


(M^  sin  a.)  q  (L1)  sin  a[sin  Q|  f  L  C^^d 


CN  ~  CNI  + 


' NC 


c 

°  ° NCO 


dx 


(6-12) 


Here,  Cdc  is  the  crossflow  drag  coefficient  for  an  infinite  circular  cylinder,  q(L')  is  a  correction  factor  for 
a  finite  length  of  L'  "calibers”,  and  d  is  the  projected  width  of  active  loading  perpendicular  tc  the  negative  of 
the  crossflow  velocity  vector.  Supersonic  leading  edge  width  on  fin  spans  are  not  included.  Side  edges  and 
subsonic  leading  edges  widths  are  included.  L'  is  a  mean  caliber  value  given  by 


L' 


(6-13) 


Cnc  Cnco  are  the  modified  Newtonian  estimates  for  the  cross  section  of  interest  and  for  a  circle  of 
diameter  d,  respectively.  Note  that  a  should  be  with  respect  to  the  Cn  =  0  angle  of  attack  for  half-plane  body 
symmetry. 
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What  is  really  being  corrected  by  the  Newtonian  ratio  is  the  different  wake.  The  impact  side  should  be 
estimated  by  Cni,  which  is  the  inviscid  normal  force  coefficient.  Hence,  the  ratio  refers  to  the  lee  side.  For  a 
body  that  is  not  symmetric  about  a  horizontal  plane,  the  crossflow  lift  for  negative  a  is  different  than  for 
positive  a. 

The  pitching  moment  is  given  by  a  similar  integration: 


C  =  C  ,  - 

m  ml 


Crfc  i\(L')sina  |sma| 


arlr 


L  (x-x')C 


NC 
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'NCO 


Here,  a  should  be  with  respect  to  the  Cm  =  0  angle  of  attack.  Here,  x'  is  the  moment  center  location,  Lr 
is  the  moment  reference  length,  and  Cm(  is  the  inviscid  moment  contribution. 

The  evaluation  of  Cnc  and  Cnco  requires  further  development.  Figure  6  illustrates  some  possibilities. 
For  half  a  circular  cylinder,  the  Newtonian  normal  force  is  1  d  Cpo.  For  the  wake  surfaces, 
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Here,  s,  is  the  arc  length  and  the  integral  is  the  over  all  wake  surfaces. 
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For  a  body  with  quarter-plane  symmetry,  Equation  (6-16)  is  0  at  0  =  n/2.  For  a  body  with  half-plane 
symmetry,  Equation  (6-16)  will  determine  the  shadow  boundary  at  0  =  0V.  For  a>0,  O<0<0V  are  the  limits  of 
integration.  For  a  <  0,  the  limits  are  0V  <  0  <  n. 


For  the  fin-body  section,  it  is  questionable  as  to  what  should  be  the  wake  surface  and  how  much  to 
integrate  over.  The  choice  taken  here  is  to  treat  the  body  and  fins  independently. 


The  value  of  d,  then,  is 


d=d.  +  2 

O 

Here  NF  is  the  number  of  fins  in  the  half  plane  and  Sn 
edge  that  is  subsonic.  Cnc/Cnco^  becomes 
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n  =  1 

=  SO-SI  for  each  fin  that  has  a  side  edge  or  leading 
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Here,  8ma  is  the  maximum  9  and  9mi  is  the  minimum  9  as  determined  by  Equation  (6-16). 

Thus,  for  a  circular  cross  section  with  a  monoplanar  wing  at  9fn  =  90  deg  and  semispan,  s,  one  obtains 
Cnc/Cnco^  =  db  + 3s  for  a  subsonic  edge,  db  is  the  value  for  a  supersonic  edge. 


In  the  application  of  Equations  (6-12)  and  (6-14),  the  values  of 
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are  only  geometry  dependent. 

In  evaluating  Cni  and  Cmi,  one  must  sometimes  reduce  the  value  of  a  to  a'  in  order  to  obtain  a 
computation  at  high  angle  of  attack.  The  control  deflection  also  must  be  reduced  for  large  control  deflections 
by  a  factor  trg.  For  configurations  without  quarter-plane  symmetry,  a  must  be  relative  to  the  Cn  =  0  angle  of 
attack. 


Pitch  plane  coefficients  are  then  adjusted  as 


CNI  =  CNI  sin  a/sina’  cosa  +  {CN  —  Cv)/f  rg 


C..,  =  C.,Tsina/sina'  cos  a  +  (Cw,  -  C...  Vt  . 
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(6-19) 
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Implementation  of  Equations  (6-19)  and  (6-20)  requires  two  separate  runs:  the  first  run  is  at  a=a'  and 
the  control  deflections  are  zero;  the  second  run  is  at  a=a'  and  the  control  deflections  are  reduced  by  trs- 
Single  prime  values  are  for  the  first  run.  Double  prime  values  are  for  the  second  run.  The  values  of  a  and  a' 
must  be  altered  for  half-plane  symmetry  as  earlier  discussed. 


The  axial  force  coefficient  prediction  requires  more  extensive  computation.  The  axial  force  coefficient 
is  assumed  to  vary  as 


CA  =  (CAO  +  CA2  Sin2  Q)  roS  Q  +  ^  5  {CADD 


tan5  +  C  .....sin  a)  cos  a 
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Here, 


{C ao +CA2sin2a)cosa  =  CAB  (o ,  8  =  0)  +  +  CA/ (a ,  8  =  0) 


(6-22) 


Four  computations  at  (a  =  0, 8  =  0),  (a  *  0, 8  =  0),  (a  =  0, 8  *  0),  and  (a  *  0, 8  *  0)  are  required  to 
approximate  the  four  coefficients  of  Equation  (6-21).  However,  the  entire  range  of  a  and  8  variation  is 
established  for  a  given  Mach  number.  The  restart  capability  limits  the  amount  of  computation  required.  Cai 
is  the  inviscid  component  of  axial  force  coefficient.  Ca  at  large  combinations  of  a  and  8  is  nonlinear  and 
difficult  to  estimate. 
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7.0  EVALU  ATION  OF  THE  NUMERICAL  METHODS 

7.1  INVISCID  POTENTIAL  THEORY 

7.1.1  Continuous  Bodies  Alone 

A  number  of  examples  of  the  computational  feasibility  of  the  implicit  Van  Dyke  model  were  given  in 
References  6  and  7.  The  significant  differences  between  first-  and  second-order  linear  potential  theory  are 
illustrated  best  by  conical  self-similar  flows. 

Figure  7  shows  the  computational  comparison  between  the  linear  potential  models  and  the  explicit 
Euler  model,  SWINT,  of  Reference  10  and  the  full  nonlinear  potential  implicit  model  of  Reference  28.  The 
reference  area  is  the  cone  base  area;  reference  length  is  the  cone  length;  and  moment  center  is  the  nose.  The 
nonlinear  potential  code  solution  was  an  unpublished  solution  provided  by  G.  Smith  of  DEI-Tech,  Newport 
News,  Va.  The  nonlinear  solutions  are  distance  asymptotic,  while  the  linear  potential  solutions  are  given  by 
an  unfactored  matrix  inversion.  The  linear  potential  solution  utilized  a  JM=  19,  KM  =  10  uniform  grid.  The 
nonlinear  solver  grids  were  finer.  Figure  7  indicates  that  the  second-order  solution  compares  fairly  well  with 
nonlinear  solutions  and  is  a  significant  improvement  over  the  first-order  solution.  Run  time  is  much  faster 
for  obtaining  a  starting  conical  solution. 

Figure  8  illustrates  another  cone  solution  for  a  body  with  a  curvature  discontinuity  in  the  cross  section. 
Stancil  refers  to  the  second-order  linear  potential  surface  singularity  code  of  Reference  29  (the  unpublished 
solution  was  provided  by  R.  Stancil).  References  are  the  same  as  before.  Note  that  the  Stancil  code  runs 
longer  than  the  nonlinear  codes  with  only  fair  accuracy.  Once  again,  the  second-order  potential  provides  a 
much  improved  solution.  The  second-order  code  utilized  a  clustered  grid  (nonuniform)  with  JM  =  21, 

KM  =  10. 

Nonconical  configurations  have  smaller  body  slopes  further  back  on  the  body. 

Figure  9  illustrates  a  computation  for  a  body  with  2/1  elliptical  similar  cross  sections.  The  body  is 
described  by 


rb(x,Q)  =  X(x)T(Q)  = 


ton (20  deg)x{  1  —0.5  x ) 
(2.5  +  1 .5  cos  28)* 


(7-1) 


Reference  area  and  lengths  are  the  same  as  before.  For  the  linear  potential  codes,  JM  =  21,  KM  =  8,  IM  =  2t 
(number  of  constant  x  step  lengths).  NCOREL  refers  to  a  full  nonlinear  potential  implicit  relaxation  code 
solver  of  Reference  30.  NANC  is  the  current  linear  potential  code.  Note  that  the  first-order  solution  is  quite 
adequate  for  this  overall  slender  configuration. 


Figure  10  shows  a  very  complex  nonsimilar  cross-section  body  for  which  a  wind  tunnel  model  with 
pressure  taps  is  described  in  Reference  31.  The  body  is  given  parametrically  as 


yJXix)  -  1  + 


1.35  T(0j)  -  1 


sin 2  n* 


sin  8, 


zb/X(x)  =  1 1  + 


1.357(0^  -  1  +  O.35/COS0J 


sin2  nx 


cosQ, 


(7-2) 

(7-3) 
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TO,)  = 


1.25  -  0.25  cos 20,  +  0.13174  (cos6,  -  cos3  0 


(7-4) 


=  land 


(7-5) 


rh  =  (y?  +  z?)* 


(7-6) 


Given  0,  one  can  solve  for  0j  by  Newton-Raphson  iteration,  utilizing  Equations  (7-2)  through  (7-5).  Cp  load 
comparisons  with  data  are  shown  in  Figure  1 1 .  References  are  as  before.  For  this  fairly  blunt  body,  second- 
order  results  significantly  improve  the  first-order  results.  The  flow  is  separated  for  x  >  0.75.  However,  the 
integrated  pressure  results  agree  quite  well  with  the  various  computations. 


The  final  continuous  body  comparison  is  shown  in  Figure  12.  The  body  (described  in  Reference  32)  has 
a  Haack -Adams  area  distribution  3/1  elliptical  cross  section.  Reference  area  is  the  maximum  area  at  x/L  = 
0.68.  Reference  length  is  the  equivalent  circular  body  diameter,  Lr  =  (4  Ar/ti)*.  Moments  are  taken  about 
x/L  =  0.6.  The  Haack-Adams  equivalent  axisymmetric  body  is  described  in  Reference  33.  Departure  from 
the  nonlinear  potential  solution  can  be  attributed  to  compressibility.  Additional  departure  from  the  data  can 
be  attributed  to  crossflow  separation.  The  NCOREL  computation  is  taken  from  Reference  30.  Body 
description  near  the  nose  was  modified  so  that  the  blunted  body  could  be  approximated  as  a  pointed  body.  For 
0  <  x/L  <  0.3,  a  third-order  polynomial  for  the  equivalent  circular  body  radius  was  used  as  described  in 
Reference  7. 


7-1.2  Discontinuous  Bodies  Alone 

Some  of  the  results  to  be  shown  here  were  described  in  Reference  8. 

Figure  13  depicts  a  comparison  with  the  data  of  Reference  16  and  an  Euler  computation.  The  Euler 
computation  starts  with  a  blunt  body  solution  computed  by  a  separate  time  asymptotic  Euler  code  followed  by 
a  computation  to  account  for  the  angle  of  attack  and  then  the  explicit  marching  SWINT  computation.  The 
two  different  Van  Dyke  second-order  potential  axial  source  and  doublet  distribution  computations  of 
References  3  and  6  produce  results  superior  to  that  shown  for  the  second-order  NANC  solution.  However, 
when  the  general  approaches  of  References  3  and  6  were  used  to  create  a  starting  plane  of  data  for  an  explicit 
marching  solution,  the  results  were  inferior  to  that  shown  for  the  current  NANC  computation.  The  nose-cone 
junction  is  a  curvature  discontinuity  that  leads  to  a  jump  in  velocity  gradient  and  pressure  gradient.  It  is 
seen  that  the  axial  force  result  is  adequate,  but  that  the  normal  force  and  pitching  moment  results  are 
degraded.  The  results  are  very  sensitive  to  the  input  gridding.  An  optimum  selection  of  a  (i.e.,  the  match 
point)  and  DO  that  controls  the  step  size  may  require  several  trial  runs  up  to  the  nose-body  junction. 

Figure  14  depicts  a  demonstration  of  the  capability  of  local  2-D  jump  relations  to  predict  the  pressure  just 
downstream  from  a  planar  discontinuity.  The  Euler  solution  was  taken  from  the  tables  of  Reference  34.  The 
•first-order  solution  deviates  significantly  from  the  second-order  solution  before  the  discontinuity  as  well  as 
after.  The  second-order  linear  jump  prediction  is  comparable  to  the  Van  Dyke  hybrid  solution  of  Reference  5. 

Figure  15  depicts  the  local  pressure  distribution  just  downstream  from  an  inlet  discontinuity.  The 
inlet  lip  lies  on  a  15-deg  ray  from  the  origin.  Compressive  turning  angles  varied  from  10.8  to  15.7  deg. 

Accuracy  deteriorates  as  oblique  shock  detachment  angles  are  approached.  Second-order  accuracy  is  still 
quite  adequate  when  compared  with  the  first-order  solution.  Reference  area  for  dCff/dix/D)  is  the  cross- 
sectional  area  just  downstream  of  the  inlet  face. 
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The  last  local  jump  solution  case  is  taken  from  Reference  14.  M«  =  1.5  on  a  long  cylinder  with  all 
disturbance  velocities  equal  to  zero  before  a  12-deg  boattail.  The  cylinder  is  at  a  =  P'  =  0  flow  inclination.  For 
this  case,  Van  Dyke  showed  that  B  =  0  gives  a  superior  result.  A  table  of  second-order  Cp  values  for  various 
combinations  of  methods  is  given  below  for  comparison  with  the  Prandtl-Meyer  value  of  CP=  —  0.293: 


B=0  B= 1 

Linear  Jump  -0.274  -0.122 

Nonlinear  Jump  -0.301  -0.223 

Other  computations  show  that  E  -  0  is  superior  for  Hares  following  a  uniform  flow  also. 

A  last  comparison  is  for  the  body  of  Figure  16.  The  first  body  section  (0<x<  1.0)  varies  from  a  10-deg 
cone  to  a  10- by  20-deg  elliptical  cross  section.  The  semi-major  axis  varies  as  a  =  x  {tan  {10  deg)  +  [tan  (20  deg)- 
tan(10deg)]x}- The  second  section  (1.0  <  x  <  2.0)  has  an  elliptical  cross-section  conical  frustum  with  b  =  tan 
(10  deg)  x  and  a  =  tan  (20deg)x.  The  first  discontinuity  leads  to  a  0.0  to -8.8-deg  expansion  (compression  is 
positive)  over  the  quarter  plane.  The  third  section  (2.0  <  x  <  4.0)  is  a  constant  cross-section  ellipse  with  b  =  2 
tan  (15  deg)  and  a  =  2  tan  (20  deg).  The  computational  comparison  is  given  in  Figure  17.  No  suppression  of 
derivatives  is  used  for  the  first  discontinuity.  The  second  discontinuity  has  expansion  turning  angles  from  -7.5 
deg  to  -20.7  deg.  The  £  derivatives  were  suppressed  for  four  steps  after  the  second  discontinuity.  Reference 
area  is  the  base  area;  reference  length  is  (4  A^/n)1/2;  and  moments  are  taken  about  x  =  2.0. 


7. 1.3  Configurations  with  Lifting  Surfaces 

A  large  number  of  computations  indicate  that  superior  results  occur  for  the  approximate  fin  boundary 
condition,  F  =  0.0,  for  all  fin  computations.  Figure  18  indicates  some  comparisons  for  supersonic  leading  edge 
jump  conditions.  The  body  shape  is  ri,  =  tan  (15  deg)  x  (1-0.5  x).  Flat  plate  fins  begin  at  x  =  0.5.  The  span,  s, 
is  0.15.  For  X  =  50  deg,  the  compressive  turning  angles  are  too  large  for  attached  oblique  shock  solutions. 

Figure  19  depicts  a  number  of  biconvex  cruciform  fin  configurations  in  the  45-deg  roll  position. 
Comparisons  with  SWINT  computations  are  made  for  various  thickness  ratios,  t/c,  and  control  deflections,  5, 
as  given  in  Table  2.  Results  for  FIN  4  are  quite  good,  considering  the  fact  that  the  pressure  oscillations  were 
strong  for  the  subsonic  leading  edge  case. 

The  next  example  is  similar  to  the  previous  examples  of  Figures  19  and  Table  2.  The  FIN  1  fins  are 
rolled  to  the  plus  position.  The  horizontal  fins  are  shifted  back  one-half  caliber  on  the  body  so  that  the  fins  are 
staggered.  The  fins  are  flat  plates.  The  fins  are  all  deflected  10  deg  (roll  control).  The  body  length  is  6.5 
calibers.  Free-stream  conditions  and  reference  dimensions  are  the  same  as  for  the  previous  case.  Coefficient 
comparisons  with  a  SWINT  computation  are  given  below: 


cA 

Cn 

Cm 

Cy 

c„ 

c. 

NANC  (1st) 

0.32 

0.74 

-1.18 

0.050 

-0.23 

1.27 

NANC  (2nd) 

0.35 

0.79 

-1.23 

0.053 

-0.25 

1.37 

SWINT 

0.33 

0.82 

-1.41 

0.086 

-0.34 

1.28 

It  is  seen  that  the  first-  and  second-order  computations  compare  favorably.  The  small  induced  side  force  and 
moment  are  sensitive  computations. 

The  next  computational  example  is  for  the  configuration  of  Figure  20  and  References  35  and  36.  Data 
and  computations  are  for  M»=  1.5.  At  this  Mach  number,  the  wing  leading  edge  is  subsonic  and  the  tail 
leading  edges  are  barely  supersonic,  but  are  treated  computationally  as  subsonic.  Figure  21  compares 
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pitching  moment  and  normal  force  coefficients  for  two  roll  positions.  In  Figure  22,  side  moment  and  side  force 
comparisons  are  shown.  Rolling  moment  is  a  more  nonlinear  coefficient.  Ataj  =  10  deg,  4>— 50  deg,  the  code 
predicts  C(  =  -0.32;  the  data  value  is  Cf  =  -0.75  Reference  area  is  the  base  area  and  reference  length  is  the 
equivalent  circular  body  diameter. 

The  next  computational  example  is  for  a  body-wing-tail  configuration  from  Reference  32  A  body-alone 
comparison  has  been  made  earlier  in  Figure  12.  The  geometry  for  the  3/1  ellipical  cross-section  body  configur¬ 
ation  is  shown  in  Figure  23.  Figure  24  shows  a  comparison  of  pitch  plane  coefficients  References  are  the  same 
as  for  the  bodv-alone  case.  Nonlinearity  for  the  moment  occurs  for  small  a  values.  A  roll  control  computation 
was  attempted  at  M«  =  2.0, 8  =  7  04  deg,  a  =  U  deg  Cf  =  0  39  for  the  computation  and  the  data  value  is  Cf  =  0  5 
The  computation  is  smoothed.  An  unsmoothed  flat  plate  fin  computation  yields  Cf  =  0  45. 

Thus  far,  all  the  computational  examples  use  analytical  description  of  the  body  geometry  The  next 
case  taken  from  References  37  and  38  required  a  numerical  input  of  data.  The  blunt  nose  variable  elliptical 
cross-section  body  is  shown  in  Figure  25.  The  spherical  nose  cap  was  input  analytically.  Numerical  data 
were  generated  from  the  right-side  columns  of  Table  3.  The  nose  cap  has  a  circular  cross  section  and  a 
discontinuity  at  the  shoulder  The  cross  sections  of  the  body  vary  from  circular  to  a  3/1  ellipse  and  back  to 
circular.  The  fins  are  the  same  us  for  the  previous  example  The  body-alone  forebody  axial  force  (including 
surface  friction)  is  0.20;  the  experimental  value  is  0. 18  for  Mr,  =  25  References  are  the  same  as  for  the 
previous  example.  Figure  26  compares  pitch  plane  body-alone  comparisons  with  data  Note  that  the  first- 
order  theory  needs  a  significant  compressibility  correction.  For  this  complex  shape,  the  second  order  solution 
departs  from  the  data  significantly  above  a  4-deg  angle  of  attack  Clearly,  a  crossflow  correction  will  drive 
the  solution  in  the  wrong  direction  for  second-order  theory  First-order  results  will  be  improved  by  crossflow- 
corrections.  The  flow  over  the  rear  portion  of  the  body  probably  separates  longitudinally  as  well  as  in  the 
crossflow  plane.  Cn  and  Cm  for  the  body-wing-tail  configuration  are  compared  in  Figure  27.  C\  comparisons 
are  good  C,„  is  a  sensitive  computation  for  xcp  close  to  x' 

7.2  HIGH  MACH  NUMBER  PREDICTION 

Figures  28  and  29  depict  comparisons  for  the  configuration  earlier  considered  in  Figure  23  and 
described  in  Reference  32  The  BWTl  comparison  considers  only  the  fin  on  the  w  indward  side.  The  BWT2 
computation  takes  into  account  both  fins  Note  that  the  local  solution  method  prediction  for  the  body  normal 
force  is  above  the  experimental,  hence,  a  crossflow  correction  would  be  inappropriate  The  same  is  true  for 
the  entire  configuration  The  seme  conclusion  is  true  for  the  moment  computation  xcp  experimentally  is 
close  to  x'  for  the  entire  a  range;  hence,  the  computation  is  very  sensitive  Departure  from  the  data  seems  to 
occur  at  high  angles  of  attack  where  the  flow  turning  angles  used  for  the  local  wedge  shock  and  Prandtl 
Meyer  solutions  on  the  lifting  surfaces  become  large 

A  slightly  different  approach  is  taken  for  the  computation  for  the  configuration  show  n  in  Figure  30  and 
taken  from  Reference  39.  Figure  31  compares  linear  aerodynamics  extrapolated  from  computations  at  a  -  5 
deg  with  data  for  M*  -  2.86.  In  this  case,  the  normal  f  ree  computation  is  below  the  data  The  moment  data 
are  typically  more  nonlinear. 

Computations  for  higher  Mach  numbers  show  the  same  trend  as  in  Figures  28  and  29  for  general 
configurations  with  lifting  surfaces,  but  closer  to  the  data  On  the  compression  side  at  higher  angles  of 
attack,  shock  detachment  is  delayed  for  higher  Mach  numbers. 
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7.3.  COMPLETE  CONFIGURATIONS  INCLUDING  VISCOUS  ESTIMATES 

7.3.1  High  Angle-of-Attack  Normal  Force  and  Pitching  Moment 

Here  computational  examples  will  be  given  for  the  quadratic  dependencies  on  angle  of  attack  predicted 
by  Equations  (6-12),  (6-14),  and  (6-19)  through  (6-21). 

Figure  32  depicts  a  comparison  for  the  missile  model  of  Figure  25  at  M®  =  2.  The  normal  force 
computation  is  quite  good.  All  inviscid  computations  were  made  at  a'  =  5  deg.  The  data  departs  from  theory 
for  the  moment  at  small  angles  of  attack  as  had  been  shown  in  an  earlier  inviscid  computational  comparison. 
The  inviscid  model  is  the  second-order  potential  theory. 

Figure  33  shows  another  comparison  for  the  3/1  elliptic  body  missile  model  at  Moo=  3.95.  The  inviscid 
model  is  the  local  solution  high  Mach  number  estimate.  Lifting  surface  contributions  to  the  crossflow  terms 
are  neglected.  The  normal  force  computation  is  significantly  improved  over  that  for  Figures  28  and  29  when 
a'  =  5  deg  is  used.  The  center  of  pressure  is  also  not  predicted  well  for  this  Mach  number. 

Figure  34  shows  a  third  comparison  for  the  3/1  elliptical  cross-section  body  model.  This  time  the  effect 
of  a  -10-deg  control  deflection  was  considered;  a'  =  3  deg  was  used  for  this  case.  The  increment  due  to  control 
deflection  was  assumed  to  be  a  constant  and  was  computed  at  a  =  0  deg.  The  actual  control  deflection  about 
the  fin  plane  was  -7.04  deg.  At  M<d  =  1.3,  the  computation  was  considerably  more  sensitive  for  the  subsonic 
leading  edges.  Again,  the  data  indicates  considerable  nonlinearity  for  the  pitching  moment  close  to  zero 
magnitude. 

The  next  computational  comparison  is  for  the  configuration  of  Figure  35  and  Reference  40.  At 
M®  =  2.5,  all  leading  edges  are  supersonic.  The  computations  for  a  S  10  deg  were  made  directly.  For  a  >10 
deg,  a'  =  10  deg  was  utilized.  At  Moo  =  2.5,  the  comparison  for  pitching  moment  in  Figure  36  is  good,  since  the 
moment  center  is  not  close  to  the  center  of  pressure. 


7.3.2  Axial  Forebodv  and!  Total  Force  Prediction 

The  first  computational  example  for  this  section  returns  again  to  the  3/1  elliptical  cross-section  missile 
model.  Figure  37  compares  forebody  and  total  axial  force  coefficient  at  a  =  0.  The  base  pressure  model  leads 
to  overprediction  of  the  total  axial  force  coefficient.  The  forebody  axial  coefficient  is  predicted  somewhat 
better.  It  is  not  clear  from  Reference  32  what  accuracy  can  be  expected  for  the  base  pressure  correction. 

In  Reference  39,  the  base  pressure  axial  force  coefficient  is  given  as  a  function  of  angle  of  attack.  A 
table  of  comparison  for  Cab  is  given  below 


M 

Cab 

Cab  (Data) 

1.6 

0.128 

0.096 

2.16 

0.098 

0.087 

2.86 

0.067 

0.063 

The  current  model  overpredicts  Cab^  Considerably  more  data  are  needed  to  modify  the  current 
estimate.  At  M® = 2.86,  the  high  Mach  number  model  was  used.  Conversion  from  Cp  to  M  uses  the  isentropic 
relation.  This  case  has  very  large  turning  angles  from  the  end  of  the  cylindrical  afterbody  onto  the  boattail. 
The  expansion  angles  vary  from  9.5  to  16.0  deg. 
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The  last  computation  in  this  section  is  for  the  3/1  elliptical  cross-section  missile  of  earlier  examples. 
Figure  33  depicts  a  comparison  between  data  and  theory  for  the  forebody  axial  force  at  zero  angle  of  attack 
with  a  -10-deg  control  deflection.  At  M®  =  1.6,  the  increment  in  axial  force  due  to  control  deflection  had  to 
extracted  from  the  zero  thickness  case.  This  increment  was  added  to  the  undePected  fin  case  forebody  axial 
force  with  nonzero  thickness. 


8.0  CONCLUDING  REMARKS 


Computational  methods  have  been  developed  that  predict  the  aerodynamic  loading  on  noncircular 
body  configurations  for  supersonic  Mach  numbers. 

At  small  incidence  and  for  the  lower  part  of  the  Mach  number  range,  considerable  effort  was  expended 
in  developing  a  finite-difference  second-order  linear  potential  model  based  on  the  ideas  of  M.  D.  Van  Dyke. 

The  implementation  of  the  second-order  model  required  some  significant  new  research.  Both  first-  and 
second-order  solutions  are  obtained  by  implicit  and  explicit  finite-difference  marching  computations.  Other 
first-order  potential  solvers  utilize  surface  singularity  models  almost  exclusively.  Jump  discontinuities  of 
supersonic  leading  and  trailing  edges  of  lifting  surfaces  and  for  planar  body  discontinuities  are  predicted  by 
local  2-D  jump  relations  based  on  the  linear  equations  of  conservation.  At  the  subsonic  leading  and  side  edges 
of  fins,  the  jump  relations  are  not  applicable  and  the  computation  would  appear  to  be  elliptic  in  character. 

The  solutions  in  this  case  also  have  a  strong  square  root  singularity  as  well  at  the  leading  edge.  It  was 
anticipated  initially  that  only  a  first-order  surface  singularity  approach  would  remove  this  severe  obstacle  to 
successful  computation.  However,  an  approximate  starting  solution  or  jump  solution  at  the  fin  edge  followed 
by  a  marching  solution  did  indeed  lead  to  useful  normal  and  even  axial  force  prediction  for  fins  with  subsonic 
leading  edges.  However,  the  pressure  coefficients  are  strongly  oscillatory  and  the  solutions  are  somewhat 
sensitive.  For  other  configurations  with  no  significant  subsonic  edges,  the  pressure  predictions  are  quite 
good.  For  blunt  nose  bodies,  the  solution  is  poor  on  the  nose,  but  overall  aerodynamic  coefficients  are  quite 
adequate.  The  entire  body  must  lie  within  an  effective  origin  of  disturbances  free-stream  Mach  cone. 

The  shearing  transformation  between  the  body  surface  and  the  free-stream  Mach  cone,  the  explicit 
boundary  advancement  model,  the  local  2-D  jump  relations,  and  the  "thin  fin”  geometry  model  are  motivated 
by  the  work  of  A.  Wardlaw  in  developing  the  SWINT  code.  The  SWINT  code  was  also  very  useful  for 
generating  comparison  computions  where  no  experimental  data  were  available. 

The  local  solution  model,  when  adapted  properly,  predicts  quite  adequate  inviscid  loads. 
Computational  times  are  very  short.  The  simple  crossflow  and  average  skin  friction  models  predict  forebody 
axial  force  and  normal  and  pitching  moment  well  when  combined  with  the  inviscid  prediction  for  the  pitch 
plane.  A  full  6-deg-of-freedom  model  at  large  incidence  would  require  vortex  generation  and  tracking 
methods. 

The  current  computer  code  utilizes  approximately  250,000  octal  storage  locations  for  a  15-by-60  grid. 
This  storage  requirement  can  be  significantly  reduced  for  the  half-plane  symmetry  case  where  a  15-by-30 
grid  is  adequate.  Computational  times  depend  on  the  gridding  sizes  and  dimensions,  body  length,  Mach 
number,  and  Mach  number  range.  Computational  times  vary  from  0.05  to  30  s  for  the  CDC  CYBER  875. 
Considerable  flexibility  in  geometry  and  free-stream  data  input  modes  is  built  into  the  code.  A  separate 
user’s  guide  is  in  preparation  as  a  companion  document  to  this  report. 
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FIGURE  1  HALF-PLANE  BODY  GEOMETRY 
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FIGURE  2  THIN  FIN  GEOMETRY 
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FIGURE  6.  WAKE  SHADOW  GEOMETRY  FOR  NEWTONIAN  CROSSFLOW 
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*n  =  amax-  xr  = 1  °- XCQ  -  0  0 


Ca 

cN 

Cm 

SWINT 

O  39 

0.26 

-0  27 

POT  (2nd) 

O  36 

0.26 

■0  26 

POT  (1*t) 

0.26 

0.13 

-0  14 

F.P.  METHOD 

0.40 

0  26 

•  0  27 

FIGURE  7.  COMPARISON  COMPUTATIONS  FOR  A  drydx  =  tan  (20  deg) 
(1-0  5  cos  8),  CONE,  M*  =  2.0,  a  =  0  0  DEG,  |i'  =  0.0  DEG 


*r  =  *max-  xr  =  "»  xco  =  o  0 


FIGURE  8.  COMPARISON  COMPUTATIONS,  HALF  CIRCLE,  drydx  =  0.4 
(0<9<90deg);  HALF  2/1  ELLIPSE,  M.  =  2,  a  =  0.0  DEG,  P'  =  0.0  DEG 


FOREBODY  4:  DEVELOPED  CROSS  SECTIONS 


FIGURE  10  NONSIMIl.AR  NASA  FORKBOUY 
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Ar  =  Amax.  Xfl  =  1.0,  XCQ  =  o.o 


cA 

cN 

C-  m 

Cy 

Cn 

Ci 

NANC(ltt) 

C  173 

0.265 

0.048 

O  268 

•0  083 

-0.0043 

NANC  (2nd) 

O  209 

0.295 

0  041 

O  333 

-0  095 

•0.0071 

SWINT 

O  208 

0  312 

0.039 

•O  316 

-0.087 

-0.0070 

INTEGRATED 

Cp  DATA 

O  216 

O  270 

0  036 

-O  319 

-0.090 

-0.0070 

—  NANC  (2nd) 

DATA  (NASA/LRC) 

*  W  =  0°) 

Cl  (0  =  180°) 


0.<  V  0.8 

X/l\ 


FIGURE  11a.  Cp  (0,180  DEG)  AND  TOTAL  LOAD  COMPARISONS  FOR  FOREBODY  4, 

M*=  1.7,  a  =  5.0  DEG,  P'  =  5.02  DEG 


0  =  270° 


—  NANC 

DATA  (NASAJLRC) 

a  (0  =  »°l 

o  10  =  270°) 


FIGURE  1  lb.  Cp  (90,  270  DEG)  COMPARISON  FOR  FOREBODY  4, 
M*  =  1  7,  a  =  5  0  DEG,  p'  =  5.02  DEG 
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dC|\|/d(x/D)  =  0.1718 
o  2nd  ORDER,  LINEAR  JUMP 
dCN/d(x/D)  =  0.2712 
^  1st  ORDER  dCpj/d(x/D)  =  0.3215 

FIGURE  14.  Cp  AND  dCft/d  (x/D)  COMPARISONS  FOR  A  CONE-CYLINDER 

DISCONTINUITY 


FIGURE  15.  Cp  AND dCN/d  (x/D)  COMPARISONS  FOR  A  CIRCULAR  INLET 

DISCONTINUITY 
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FIGURE  19  BODY  TAIL  CONFIGURATIONS  IN  CALIBER  DIMENSIONS 


WING  LEADING  AND  TRAILING  EOGES:  DOUBLE  WEDGE 
TAIL  FIN  LEADING  EDGE  DOU3LE  WEDGE 


FIGURE  20.  2/1  ELLIPTICAL  CROSS-SECTION  BODY  CONFIGURATION 
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-  NANC 

DATA 

O  tf>  -  0  deg  ROLL 
Q  4>  =  60  dag 


TOTAL  ANGLE  OF  ATTACK  (deg) 


figure  21  cN  and  cm  comparisons  for  a  2/1  elliptical  cross 

SECTION  MISSILE,  Mx=  1.5,4>  =  0  DKG,  <$>  =  50  DEG 


— 

.  NANC 

DATA 

□ 

<t  =  90  deg  (/?  *“  0,  a  =  0) 

12  ot  (deg) 


FIGURE  22  C y  AND  Cn  COMPARISONS  FOR  A  2/1  ELLIPTICAL  CROSS 
SECTION  MISSILE,  Mx.  -  1  5,  <J>  =  90  DEG,  <J>  =  50  DEG 
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FIGURE  26.  Cn  AND  Cm  COMPARISON'S  FOR  A  VARIABLE  ELLIPTICAL 
CROSS-SECTION  BLUNTED  BODY.  M„=  2.5 
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FIGURE  27  CN  AND  Cm  FOR  A  VARIABLE  ELLIPTICAL  CROSS-SECTION 
BWT  CONFIGURATION,  M«  =  2  5 
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FIGURE  37  FOREBODY  AND  TOTAL  AXIAL  FORCE  COMPARISONS  FOR 
A  3/1  ELLIPTIC  BODY  MONOPLAN'AR  MISSILE  MODEL,  a  =  0  DEG 


FIGURE  38  FOREBODY  AXIAL  FORCE  COMPARISON  FOR  A  3/1  ELLIPTIC 
BODY  MISSILE  MODEL,  a  =  0, 8f=  -  10 
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TABLE  l  LOCAL  SOLUTION  PREDICTION  METHODS 


Cross-Sectional  Shape 

Body  Section 

(Windward  Plane) 

(Leeward  Plane) 

Circular 

Nose 

Dahlem-Buck 

Dahlem-Buck 

Body 

Inclined  Cone 

ACM  Emperical 

Noncircular 

Nose 

Dahlem-Buck 

Dahlem  Buck 

Body 

Tangent-Cone 

ACM  Emperical 

Fins/Wing 

Tangent-Wedge 

Prandl-Meyer 
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irm  »r*  vw  wr «  4  m*  *-  «mum  -M  vw«f.XMR  IT*  Vf  1^ 

i 


TABLE  2.  COMPARISON  COMPUTATIONS  FOR 
BODY/FIN  I  CONFIGURATIONS,  M„  =  2.0,  a  =  5  DEG, 
x'  =  3.0  CALIBERS  FROM  NOSE,  Ar  =  CYLINDER  CROSS- 
SECTION,  Lr- CYLINDER  DIAMETER 


FIN  1 


Case 

CA 

CN- 

r 

0,  t/c  =  0.0 

NANC  (first) 

0.091 

0.81 

-1  09 

NANC  (second) 

0.097 

0.83 

1  09 

SW1NT 

0.096 

0.83 

1  13 

=  0,  t/c  =  0.1 


FIN  1 


NANC  (first) 

0  21 

0.78 

-1  03 

NANC  (second) 

021 

0  77 

-0  90 

SWINT 

0  19 

0.77 

-0.96 

-10.0,  t/c  =  0.0 

NANC  (first) 

0  21 

-0.26 

1  52 

NANC  (second) 

0.21 

-0.24 

1  54 

SWINT 

0.19 

-0.26 

1.54 

-10.0,  t/c  =  0.1 

NANC  (first) 

0.30 

-0  22 

1.45 

NANC  (second) 

0.30 

0.17 

1.33 

SWINT 

0.29 

-0.18 

1  29 

FIN  2 


=  0,  t/c  =  0.1 


NANC  (first) 

0  23 

1.01 

NANC  (second) 

0.23 

1  00 

SWINT 

0.21 

1.00 

FIN  3 


=  0,  t/c  =  0.1 


NANC  (first) 

0  24 

0.98 

1  65 

NANC  (second) 

0.25 

0.97 

-1.53 

SWINT 

0  21 

0.88 

-1  39 

FIN  3 


=  0,t/c  =  0.0 


FIN  4 


NANC  (first) 
NANC  (second) 
SWINT 

0.091 

0.097 

0.096 

1.02 

1.03 

0  96 

1  79 
-1.77 
-1.65 

0,  t/c  =  0.0 

NANC  (first) 

0.091 

1.04 

-1.18 

NANC  (second) 

0.097 

1.18 

-1  37 

SWINT 

0.096 

1.13 

-1  50 
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TABLE  3.  SEMI-MAJOR,  SEMI-MINOR  AXIS  VARIATIONS 
FOR  VARIABLE  ELLIPTICAL  CROSS-SECTION  BODIES 


Sharp-Nose  Body 


Blunt-Nose  Bodv 


x/L 


a/L 


b/L 


x/L 


a/L 


b/L 


0  0000 


0.0000 


0  0000 


0.0000 


0.0000 


0.0000 


0.0100 


0.0071 


0.0024 


0  0216 


0.0217 


0.0217 


0.0200 


0.0118 


0  0039 


0  0258 


0.0222 


0.0219 


0.0250 


0.0140 


0  0047 


0.0799 


0.0337 


0.0257 


0.0500 


0  0234 


0.0078 


0.1340 


0.0452 


0.0287 


!  0  0750 


0.0314 


0.0105 


0.1881 


0.0564 


0.0311 


0.1000 


0.0388 


0.0129 


0  2423 


0.0673 


0.0332 


0.1250 


0  0455 


0.0152 


0.2964 


0.0777 


0.0349 


0.1500 


0.0518 


0.0173 


0.3505 


0.0875 


0  0365 


0.2000 


0.0633 


0.0211 


0.4046 


0.0965 


0.0378 


0.2500 


0.0737 


0.0246 


0  4587 


0.1048 


0.0389 


0  3000 


0.0831 


0.027/ 


0.5669 


0.1177 


0.0405 


0.3500 


0.0915 


0.0305 


0.6210 


0.1220 


0.0410 


C.4000 


0.0991 


0.0330 


0  6481 


0.1233 


0.0412 


!  0.4500 


0.1059 


0.0353 


0.6800 


0.1237 


0.0412 


0.5000 


0.1118 


0.0373 


0  7000 


0.1229 


0.04 '3 


0.5750 


0  1188 


0.0396 


0.7292 


0.1190 


0.0418 


'  0.6000 


0.1206 


0.0402 


0  7500 


0.1 166 


0.0421 


0.6250 


0.1221 


0.0407 


0.7833 


0.1093 


0.0431 


0  6500 


0.1231 


0.0410 


0.8000 


0.1066 


0.0435 


0.6800 


0  1237 


0.0412 


0.8314 


0.0965 


0.0452 


0.7000 


0.1234 


0.0411 


0.8500 


0.0945 


0.0458 


0.7292 


0.1221 


0.0407 


0.8716 


0.0867 


0.0476 


0.7500 


0.1213 


0.0404 


0.9000 


0.0813 


0.0492 


0.7833 


0.1188 


0.0396 


0  9250 


0.0748 


0.0514 


0.8000 


0.1180 


0.0393 


0  9500 


0.0687 


0.0540 


0.8314 


0.1146 


0.0382 


0  9750 


0.0635 


0.0567 


0.8500 


0.1139 


0.0380 


1.0000  I  0  0595 


0.0595 
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A 

ar 

B  =  0,1 

C 

Ca 

Cab 

Caf 

Cf 

Cf 

Cm 

Cm 

Cn 

CN 

cp 

CY 

D 

DO 

E 


NOMENCLATURE 

Free-stream  speed  of  sound 

ul  (i  qi2  —  (2-y)/6  U|2  -  1]  +  coso  cosp  [q;2  +  qj 

qi  +  (y-l)/2  M*2  Ui),  second-order  gradient  component. 

Reference  area 

Body  boundary  condition  factor 

Source  vector  in  conservation  equation  (see  E) 

Axial  force  coefficient 

Base  axial  force  coefficient 

Skin  friction  axial  force  coefficient 

F/(Ak  Qp),  general  force  coefficient 

Roll  moment  coefficient 

Pitching  moment  coefficient 

M/(Ar  Lr  Qq),  general  moment  coefficient 

Yawing  moment  coefficient 

Normal  force  coefficient 

(p/p®-  *)/(y/2  M»2),  pressure  coefficient 

Side  force  coefficient 

Circular  body  diameter  or  equivalent  circular  body 
diameter 

Af,  at  the  body  boundary 

Conservation  vector  (two  or  three  components). 
Unified  almost  conservation  form  vector  equation  for 
flow  field,  body  surface,  and  fin  surface  is 
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» 


Ma 

n 


Free-stream  Mach  number 

Normal  distance  from  a  body  cross-section  boundary 
curve  (in  the  plane  of  section)  or  unit  normal  from  a 
solid  surface 

Local  pressure 
Free-stream  pressure 


All  q  values  below  are  nondimensional 


q 

Oc 

% 

<C 

qe 

% 

% 

% 

% 

<c 

<u 

y 

Qd 

r 

Tb 

f|,  r0 

R 


Perturbation  velocity  vector 

Crossflow  component  of  perturbation  velocity  vector 
M.*4  (qi(2)(2-y)A2p2) 

v\  qy^/2  +- t'x  (q^  qi  +  qi2)  +  (y-1)  M«2/p2  [vJ2  (t>i2  - 

W l2)  -I-  W x  L  ]  Wl  1 

wi  qj 2/2  +  lu m  (q^  •  q>  +  qi2)  +■  (y-1)  M*2/p2  l  -  wJ2 
(i’l2  -  w  12)  +  v ^  cq  tail 
First-order  perturbation  vector 
First-order  crossflow  perturbation  vector 
Second-order  perturbation  vector 
Second-order  crossflow  perturbation  vector 
Free-stream  velocity  vector 

Crossflow  component  of  free  stream  velocity  vector 

Total  velocity  vector 

Px  y / '  Mr2,  dynamic  pressure 

General  cylindrical  coordinate  radius 

Body  radius 

Minimum  and  maximum  rudiul  distance  on  a  fin  for  a 
vulue  of  x 

r  lx 
r^/x 


MS 
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RMI 

RMFN 

RNF 

RNl, 

s 

SMA 

S 

Sw 

t 

t* 

T 

Ui,u2 

u2e 

LUx 

u 

^l.L'2 

l'2e 

v» 

V 

ysl  i  vs2 
W  i,  1^2 

w2v 

wr 

V/ 

x 

x'  --  XCG 


Minimum  fin  radial  distance 

Reynolds  number  per  Mach  number  per  foot  per  million 
• —  px  a  ®/(p®  106) 

Reynolds  number  per  foot  =  RMFN  (106)(MX) 

Element  Reynolds  number 
r  -  RMI  or  curve  arc  length 
Maximum  s  for  a  fin 
See  E 

Wetted  surface  area 
Fin  thickness 
Fin  slope 
See  E 

x  components  of  q5  and  q2 
u2  +  M<*,2/p2  a 
cosa  cosp,  x  component  of  q™ 
x  component  ofQ 
r  component  ofq\,  q2 
t’2  -  Mx2  qr 

sinfsinO  +  sina  sinP  cos8,  r  component  of  qx 
r  component  cf  Q 

First  and  second-order  tangential  velocity  components 
0  components  of  qi ,  q2 
w2  —  Mx2  qa 

sinP  cosO  sina  cosp  sint),  0  component  of  q.,. 

0  component  of  Q 
Distance  from  nose 
Moment  center  from  the  nose 
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XMI 

XMA 

^  * 

'■  R 

y 

z 

a 

P 

P' 

V 
8f 
< 

e 

0r 

K 

A 

* 

P„ 

o 

<t> 

V 

V 
Vc 


Minimum  x  for  a  fin 
Maximum  x  for  a  fin 
See  Lr 

Cartesian  distance  to  right  looking  downstream 
Cartesian  distance;  up  looking  downstream 
Angle  of  attack 

(M„2  m 

Sideslip  angle 
Heat  capacity  ratio 
Fin  deflection 

Clustering  transformation  for  £  variable 

Angle  distance  from  z  axis 

Angie  location  of  a  fin 

Curve  radius  of  curvature 

Fin  sweepback  angle 

Free-stream  viscosity 

Shearing  transformation 

Free-stream  density 

Angle  from  fin  midplane  to  fin  surface 

Angle  clustering  transformation  variable,  roll  angle 

9  -  Of  -  o 

Gradient  operator 
Crossflow  gradient  operator 
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